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PREDICTION OF UNSTEADY AERODYNAMIC LOADINGS CAUSED 
BY LEADING EDGE AND TRAILING EDGE CONTROL SURFACE 

MOTIONS IN SUBSONIC COMPRESSIBLE FLOW COMPUTER 

PROGRAM DESCRIPTION 


By 


M. C. REDMAN AND W. S. ROWE 


SUMMARY 


A digital computer program has been developed to calculate 
unsteady loadings caused by motions of lifting surfaces with 
leading edge or trailing edge controls based on the subsonic 
kernel function approach. 

The pressure singularities at hinge line and side edges have 
been extracted analytically as a preliminary step to solving 
the integral equation by collocation. 

The program calculates generalized aerodynamic forces for user 
supplied deflection modes. Optional intermediate output includes 
pressure at an array of points, and sectional generalized forces. 
From one to six controls on the half span can be accommodated. 


Page Intentionally Left Blank 



1 . 0 INTRODUCTION 

This document describes in detail the design and usage of the 
FORTRAN IV digital computer program, RHO IV. The RHO IV program 
was written as an engineering tool to be used in calculation 
of unsteady aerodynamic loadings on lifting surfaces ; with leading , 
edge or trailing edge control surfaces in compressible subsonic 
flow, per the analysis techniques presented in reference 1. 

Features of this program include: 

• Modal input in the form of surface deflections at 
arbitrary points - Input from disk file, tape, or cards 

• Calculation of unsteady pressures at arbitrary points 
on the lifting surface planform 

• Calculation of sectional generalized forces at arbitrary 
spanwise locations on the lifting surface planform 

• Calculation of generalized forces 

• Optional saving of unsteady pressures, sectional or 
total generalized forces for subsequent analysis 

• Optimization of computer time through the capability 
- to save and later reuse pressure/kernel influence 

coefficient matrices (C-matrices) 

• Capability to provide for airfoil thickness corrections 
using velocity profile modifications supplied by the . 
user 

• Capability to analyze up to six separate closed gap 
control surfaces 

• Capability to analyze coupled main surface and control 
surface modes 

Included in this document are: 

• Description of equations used in the program 

• . Description cf variable length storage requirements 

• Description of user 1/0 and scratch file formats 

• List of program limitations 

• Description of computer program usage 
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Data stacking procedures 

Description of program output - normal and diagnostic 
Sample problem input/output 

Description cf program structure and routines used 


2.0 DISCUSSION 

2.1 GENERAL REMARKS 

This section will present a description of the nomenclature and 
analysis as it appears in the RHO IV program. Although some 
information is given with each section, no attempt is made to 
develop or reference the sources of development of the theory 
involved. A full discussion of the'latter is given in reference 
1 of this document. 

As described in reference 1, the problem of identifying the 
unsteady aerodynamic loading on a lifting surface, without downwash 
discontinuities may be written as a boundary value problem. 

w(x,y) = l //AP (£,n)*K(x-£,y-n,k,M)d£dn (2.1-1) 

4mp\F 


In equation 2.1-1, the left hand side, W(x,y) , is the K inematic 
D ownw ash, or effective angle of attack of the surface, due to 
the structural vibration mode. The right hand side of the 
equation is the mathematical downwash which is derived from the 
surface integration of unsteady pressure times an aerodynamic 
influence function. The latter, called the Kerne l Functi on, 
is dependent upon geometry, reduced frequency, and Mach number. 
The unsteady pressure is unknown; however, knowing the physical 
characteristics of loading, the unsteady pressure may be 
approximated by a linear combination of A ssumed P ressure Term s 
which will satisfy loading characteristics. ^ 

AP r (C,n) = 4pV 2 /S‘-n 2 lajApj (£ ,ri) , j=l,m (2.1-2) 

If the Downw a sh Integral Sauation, 2.1-1, is written for a number 
of descrete points on the surface, (D ownw ash Points , or 
collocation points, or control points) , the resulting Co mplex 
L inea r System of E qua t ions may be expressed in matrix form as, 

{W (x,y) } = [C] {a} 


where the elements of the C -matri x are 


(2.1-3) 



C. . = 1 /7Ap . (£,n) •K(x i -C/y i -nfk,M)dCdn 

■* IT - '• 


(2.1-4) 


(Note that the C-matrix terms are independent of structural 
mode.) For a simple lifting surface problem, solution of 2.1- 
3 for the unknown coefficients of the assumed pressure terms, 

{a^}, allows one to calculate the unsteady pressure at any point 
on J the surface or integrate the unsteady pressure times modal 
displacements to give generalized force. 

With the introduction of control surface motion relative to the 
remainder of the lifting surface (main surface) the problem 
becomes somewhat more difficult. In particular, the kinematic 
downwash distribution or sheet will contain a step discontinuity 
at the control surface with respect to the main surface. The 
use of additional assumed pressure terms with unknown coefficients 
to match the discontinuous boundary condition is prone to numeric 
difficulties. However, a pressure term associated with control 
surface rotation may be written which will give the required 
jump at the hinge and control surface side edge and which has 
a known coefficient. Thus 


W* (x,y) . =• 


4irpV 


ff AP_ (£,n) *K (x-£,y-n,k,M) d£dn (2. 1-5) 

T s 


where W* (x,y) will have the same jump discontinuity at the hinge 
and side edge as W(x,y) and will be relatively smooth away from 
the hinge. The control surface unsteady pressure term in equation 
2.1-5 is a function of the control surface rotation at n, 9(h), 
i.e„, is dependent upon the structural mode. The control surface 
rotation is approximated by a cubic equation. 


0(n) = Ao + +. A iR cs 


( 2 . 1 - 6 ) 


where = (n-y i > / (Y^Y^ 


y^ = Inboard side edge of control surface 
y Q = Outboard side edge of control surface 
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This representation of control surface rotation should suffice 
for a broad range of control surface twist. The expression for 
the assumed control surface pressure term is then 


ap (£ / n) = 4pv 2 /s^- n z za. . A P . ( c / n ) 

s 1 J J 


(2.1-7) 


If equation 2.1-5 is written for the downwash points and expressed 
in matrix form. 


{W* (x,y) } = [C*]{A} 


( 2 . 1 - 8 ) 


where the elements of the Control Surface c-matrix are 


CJj = JL //Apj (?,n) •K(xv-C,y i -h/k,M)d5dt) 


(2.1-9) 


and significantly, because of the polynomial representation of 
control surface rotation, are not dependent upon structural mode. 
If the kinematic downwash is modified by removing any 
discontinuity due to control surface rotation. 


W(x,y) = W(x,y) - W* (x,y) 


the resulting residual dow T nwash sheet, W(x,y) , which is smooth, 
may be used to solve for the unknown coefficients of the assumed 
main surface pressure terms, {a-}. The total unsteady pressure 
is then a combination of main surface pressure and control surface 
pressure. 

AP(£,n) = AP{S,n) + ap (£,n) (2.l-/i) 

Ju S 

Thus unsteady pressure may be calculated at any point on the 
surface or integrated to produce generalized forces as in the 
simple lifting surface problem. 



2.2 NOMENCLATURE 

The RHOIV program works with dimensional coordinates (£,n) and 
non-dimensional coordinates (£,tv) . The b Q reference length used 
in reference 1 for k value and non-dimensionalizing of all 
coordinates is used in RHOIV cnly to arrive at k = w/v from the 
user input k = b 0 w/v. Note that differences in nomenclature 
between reference 1 and this section are noted parenthetically 
with the symbols. 

SYMBOLS DESCRIPTION . ■ 


A(n;m, j) 
a (m, j ) 


b 0 

bfl) 



(x,y;m) 


C* (x r y,n;m) 


C 

Cx#C 2 


E i » e 2 


Cubic coefficients (m=T,...,4) in n for 
mode j , for control surface n. 

Undetermined coefficients of assumed main 
surface pressure terms for mode j. 
(m=1,... # no. of downwash points) 

Reduced frequency reference length. 

Local planform semi- chord, 
b (n) = o.5[£ t (n) - ?! (n) ] 

C-matrix terms for downwash point (x,y) 
associated with assumed main surface 
pressure terms. 

(m=1,.. . ,no. downwash points) 

C-matrix terms for downwash point (x,y) 
associated with control surface n 
pressure terms. 

(m=1 ,4) 

Chordwise integral of g (£ r n) K (x 0 , y 0 ,k,M) , 
Eqn. 2. 3-2. 

Coefficients associated with control 
pressure expressions, 

Eqns. 2.5-11, 12. 

Chordwise modification functions 
associated with control pressure 
expressions, Eqns. 2.5-16, 17. 


f 



(n; i) 


"Spanwise" portions of a pressure term. 
(i=1,...,no. of chords, main surface 
analysis) (i=1,...,4, control surface 
analysis) 
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SYMBOL 


DESCRIPTION 


F (x,y,n) 
G(x>y».n) 


' G „-, » G C 7 » G T l> G rV G 


'IS' Si' S2 

g(£,n;j) 

H 


Ll * L2 ' AT 


Ii 

j 

K(x 0 »yo»k,M) 

K ns (x °' yo ' k ' M) 

K s <Xa,y 0 ,k,M) 

*x 

k 

* r [k] 

Lx 

L i / L 2 

M 


M(5-x s# n-y s ) 


Portion of dcwnwash integral expression, 
Eqn. 2. 3-6 . 

Portion of dcwnwash integral expression, 
Eqn. 2.3-3. 

Portions of pressure expression 
associated with full chord control. 

"Chordwise" portion of pressure terra, 

(j= 1 , . . - , no. pts. /chord, main surface) 
(j=1, control surface analysis) 

Spanwise modification function 
associated with control pressure 
expressions, Eqn. 2.5-11. 

Integration limit in kernel function 
evaluation. 

or i-th displacement node, or i-th 
"spanwise' 1 pressure term. 

Modified Bessel function. 

j-th displacement mode, or j-th 
"chordwise" pressure term. 

Full kernel expression. 

Non-singular portion of K. 

Singular portion of K. 

Modified Bessel function 

Reduced frequency, k = u)/V. 

Reference reduced frequency, k =b 0 to/v. 

Struve function. 

Portions of pressure expression 
associated with full chord control. 

Mach number 

Portion of pressure expression 
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SYMBOL 


m 


N (C-x s# rry s ) 


AP(5,n; j) 


AP r U,n;j) 


AP_ (C»n#n; j) 

CAP ae (5 f n)l 


Apcc#n;ni) 


Ap(£,n#n;m) 


q 

Qi 

Q (i# j) 
Q r (i* j) 

Q (n;i, j) 

s 

$ r (i , rc) 



DESCRIPTION 


associated with partial chord control. 

A pressure term number. 

Portion of pressure expression 
associated with partial chord control. 

A control surface number. 

Total change in pressure for mode j at 
point (£, n) • P (lower) -P (upper) 

That portion of AP associated with the 
regular (assumed main surface) pressure 
terms. 


That portion of AP associated with the 
pressure terms for control surface n. 

The m-th assumed main surface pressure 
term value at (£,n) • 

The m-th pressure term associated with 
control surface, n, value at (£,n)* 

Dynamic pressure, q = 0. 5pV 2 . 

Portion of pressure expression 
associated with full chord control. 

Generalized force (generalized unsteady 
aerodynamic coefficient) for displacement 
mode i, pressure mode j. 

That portion of Q associated with the 
regular (assumed main surface) 
pressure terms. 

That portion of Q associated with the 
pressure terms for control surface n. 

Surface integral of Z (£,n;i) AP * 
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SYMBOL 


DESCRIPTION 


Q (n;i#j) 


Q In sir j) 


Sectional generalized force (sectional 
generalized unsteady aerodynamic 
coefficient) for station n# displacement 
in mode i , pressure in mode j. 

That portion of Q s associated with the 
regular (assumed main surface) pressure 
term. 


Q (n#n;i,j) 


That portion of Q 5 associated with the 
pressure terms for control surface n. 

Chordwise integral of Z(£,n;i) Ap(£,ri;m) . 

A modified distance between points in 
kernel and control pressure expressions. 

Planform semi span 

Time 


u [V t ] 


W(x,y; j) 

W*{x,y,n; j) 


W(x # y; j) 


Local streamwise velocity. 

Remote freestream velocity. 

Kinematic downwash at (x r y) for mode j. 

Mathematical downwash at (x,y) for 
mode j due to control surface n 
pressure terms. 

Residual downwash at (x,y) in mode j ; i.e. 
W with all control surface 
discontinuities removed. 


x [b 0 x] 


X ,x . , x 
S X o 


Local normal (z) velocity. 

Downwash point chordwise coordinate 

Control surface hinge corner 
(i = inboard) , (o = outboard) 


y [b 0 y] 


Downwash point spanwise coordinate. 
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SYMBOL 


y #y. »y_ 

s i o 


z (x,y,t) [Z ] 


Z (x,y;i) 


6 

% 

3l 

n [ b on] 


n 




Q(n) [0 H ] 


5 £bo€] 

ICn) 

C c (n) 

?! in) 



^t (n ) 


DES CRIPT I ON 

Control surface side edge (i=inboard), 
(o=outboard) . 

Time dependent modal displacement 
normal to surface. 

Modal displacement normal to surface 
at (x,y) in mode i. 

/T=rr* 

/^ z +tan z A H 

/$^+tan 2 A L 

Dimensional spanwise coordinate. 

Non-dimensional spanwise coordinate, 

_n = n/s. 

Non-dimensional spanwise coordinate 
referenced to control surface span, 

ics * (n-Yi>'(y 0 -yi>- 

"Chordwise" non-dimensional pressure 
term coordinate, 0 = cos“l (-£) 

Control surface streamwise hinge 
rotation. 

Control surface hinge sweep. 

Control surface leading edge sweep. 
Dimensional chordwise coordinate. 
Non-dimensional chordwise coordinate, 

i *■ EC-e m «n) ]/b fn) - 

Control surface hinge value. 

Planform leading edge value. 

Planform midchord value, 

= C? t (n)+C l {n) ]/2 
Planform trailing edge value. 
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SYMBOL 


DESCRIPTION 


Fluid mass density. 

"Spanwise ,, non-dimensional pressure 
term coordinate, <j> = cos _i (ri) 

Circular frequency of oscillation. 
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midchord 



1 


2.3 DOWNWASH INTEGRAL EXPRESSION 

The downwash integral expression, which relates kinematic and 
mathematical downwash, is 

S . . ... ' 

w (x,y) =1 //s 2 -n 2 [f (n)c + G(x,y,y)/y* - .G ' (x,y ryl/y^J-dn 2.3-1 

TT -S 

+ G (x,y ,y) + yC (x,y,y) 

where C is the integral of the product of the chordwise pressure 
term and the kernel function, 

Xl. (n) (n) 

C = / g(C,h)K (x ,y n ,k,M)dC + / g (C/h)K (x n ,y f k,M)dC 2.3-2 
(n) ns 0 0 qcn) s 0 0 

and G(x,y,y) and G' (x,y,y) are related to the evaluation of the 
dipole singularity; see reference 1. 

G(x,y,n) = f(n)F(x,y,n) 2.3-3 

G*(x,y,n) = 3G(x,y,n) /3h = f Mri)F{x„y,n) + f (n) F’ (X,y,n) 2.3-4 

* G (x,y,y) = Lim G(x,y,g) 2.3-5 

n-*r 

where Ap(£,n) = f (rj) g (C,h) is the loading function, and 


F(x,y,n) 

= / g(?/n) [1 + x n /Vx* + 6 i yn]e ixx 0d5 

2.3-6 


^(n) ° 


F (x,y,y) 

x ikx 

= 2-/g(C,y)e 1KX 0dC 

2.3-7 


h (y) 


F* (x,y,n) 

= 9F(x,y,n) /3n 

2.3-8 

F' (x # y,y) 

X ikx 

= 2*/{3g (?,y) /3n] e 1KX 0d^ 

2.3-9 


?l (y) 

Note that integration by parts, or a similar approach, is required 
in equation 2.3-9 for those terms, 3g(£,y)/9ri, which contain 
a singularity in the interval, see reference 1. 


Note also, that the spanwise integrand associated with the 
singular kernel, i.e. , 


S t (n) 

✓s 2 -n z [f (n) / g (C'»n)K (x n ,y n ,k,M)d£ + G(x,y,y)/y^ - 
C x (n) s 0 0 0 


c (x,y,y)/y 0 ] 


is of the form, 

ln|y 0 /S | (regular function) + (regular function) 

which requires log quadrature to be used around a downwash chord 
for that portion of the spanwise integrand. 



2.4 KERNEL FUNCTION 

The kernel function is an aerodynamic influence function relating 
the induced normal velocity at an arbitrary field point to a 
unit loading on a surface at some other point. In the case of 
a flat plate lifting surface, only the planar portion of the 
kernel is used. This may be written, 

K s (xo,y 0 ,k,M) = k 2 e‘ lkxo {-(R+kxo)/(Rk 2 y§) + i/R - ( kx 0 -MR ) / ( 2 e 2 R ) 

-In j (R-kxo)/(2-2M) j/2| 

where R = /k 2 x 2 + g 2 k 2 y 2 


The expression 2.4-1 has been shown to contain a number of 
singularities which cause numerical integration to be extremely 
expensive. The singularities have been identified and may be 
analytically subtracted from the full expression yielding a 
nonsingular function. The form of the singular portion of the 
kernel is 


K(x 0 ,y 0 ,k,M) = -e lkXo l /[e^V/T^^dA 

y o 3y <T~ 


2.4-2 


where h = [x () -M» / x 2 + yjp/B : 


which reduces to eqn. 2.4-3 when k=0. 


M x o’ y o*°’ M ) = -{1 + x n //x J + 


o ' ^ o 


2.4-3 


The singular portion is integrated separately from the nonsingular 
portion. Because of its singular nature, it requires a large 
number of evaluations when numerical integration is being 
performed; however, it is relatively inexpensive to evaluate.. 

The nonsingular function is slightly more expensive to evaluate 
than the full kernel; however, because of its regular nature, 
it is evaluated much less during numerical integration. 

Two forms of the nonsingular kernel are used. The first, Watkin's 
formulation, is faster to calculate (and numerically sufficiently 
accurate) for values of kJy o |>1.0. The second, Rosel's 
formulation, requires longer to calculate (particularly as k|y Q | 
becomes large), but is numerically more accurate for k|y o |<1.0. 


For Watkin’s form, defining 





2.4-4 


The form of the full kernel expression is 

K(x 0 ,yo,k.M) * k 2 e ' kx 0 { - K ! ( k | y 0 1 ) / k I y o [ - i .5»[r, (k|y 0 | }-L \ (k |y 0 1 )]/k|y 0 | 2 _ 4 . 5 

+i/k|y 0 l- kx 0 e ,kh /fc 2 y8 ■ C 

The two expressions, Eqns. 2.4-2 , 5 are combined to give the 
ncnsingular form. 

K a (Xo»yo#k#M) = K{x 0 ,y 0 ,k,M) - K (x 0 ,y 0 ,k,M) 2.4-6 

ns .O’ 

In Eqn. 2.4-4, the term x/ / 1 + t 2 is approximated by an exponential 
series, reference 6, which may be integrated analytically. 

For Rosel's form, the nonsingular kernel is written directly. 


K (x 0 , y 0 ,k,M) = k^e" 1 kx ° { /[(e 1 X -1 -i x+x 2 /2 )/ (x 2 +k 2 y§ ) 3/2 ]dx 2.4-7 



where the singular terms {the same as in Eqn. 2.4-2} are already 
removed. The exponentials in the integrals are written as 
infinite series and integrated analytically. The resulting 
infinite series of terms, which may be calculated in a recursive 
manner, are truncated when a predetermined conversion criteria 


is met. 
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2.5 LOADING FUNCTIONS 

The loading functions used within RHOIV are of two types. The 
regular, or main surface, pressure terms are used to match the 
regular boundary condition associated with a simple lifting 
surface or the residual boundary condition associated with a 
lifting surface with controls when the control discontinuities 
have been removed. The singular, or control surface, pressure 
terms are used to match the c hanc e in boundary condition at a 
control surface hing e or side edge. 

The regular pressure is assumed to be of the form 

AP r (£,n) = 4pV 2 /S 4 -n i ZAp (£,n;m)a(m) 2.5-1 

TCI 

= 4pV 2 /s’ T -ri" T {Ap (C,n;m) } ’ {a(m) } 

where a(m) are unknown coefficients. The term /S 2 -n 2 has been 
included in the coefficient to simplify evaluation of the dipole 
singular portion of the downwash integral expression. 

The assumed main surface pressure terms, Ap (£, g) , are themselves 
composed of a "spanwise" and a "chord wise” term, e. g. 

Ap(£,n;m) = f (g;i)g(C,n; j) 2.5-2 

The set Ap _ ( C , n ) “is composed of a 11 combinations of spanwise, 
f (p) , and chordwise, g (£, n) , terms. For the determined case, 
the number of spanwise terms is equal to the number. of downwash 
chords (=NSPT) , and the number of chordwise terms is equal to 
the number of points per downwash chord (=NCPT) . 

The spanwise terms used are 

f(n;i) = sin[ (2»i-N) <p]/Y5 7r= r\ T r i = 1,-..,NSPT 2.5-3 

where <j> = cos- 1 (n) , n = n /s 

N = 1 for a symmetric analysis, 

=0 for an antisymmetric analysis 

The chordwise terms used are 

g(£,n;i) = cot (6/2) 

g(£#n:j) = sin[ { j- 1) 9], j = 2,3 ,...,ncpt 
where 9 = cos -U-£) , £ = (£ - £ m (n) ) /b (n) 


2.5-4 


Note that all Ap(£,n) go to zero at the planform trailing edge 
and tip as the square root of the distance., The terms which 
are not associated with g(£,n;1) go to zero at the leading edge 
in a similar manner. The terms which are associated with g(£,n;T) 
have the inverse square root singular form which is required 
at the leading edge. Note also that f (0; i) is plus or minus 
one for a symmetric analysis, and zero for an antisymmetric 
analysis. 

The (i, j) combinations of Eqn. 2.5-2 are ordered 
{ (i=1 , . . , ,NSPT) , j=1,. .NCPT) . 

The singular pressure expression associated with a control surface 
is of the form. 


AP C U,n) = pv 2 0(n)g(S,n)A 


2.5-5 


where as indicated in Eqn. 2.1-6 the streamwise control rotation, 
0 (q) , is represented as a cubic. 


0(n) = £A(m)iT™ 1 

i — cs 

^ = (n-yi)/y 0 - yi ) 

Defining the “spanwise" pressure terms to be. 


2.5-6 


‘ f(g;i) = / 4 tt/S 2 -ri z 

“C O 

Equation 2. 5-5 becomes. 


2.5-7 


AP(£,ti) = 4 P V2 ✓Ssr- T1 ^{AF(C.n;n')} ' (A(m)} 

S' 

where Ap (£,n;m) = f(n;i)g(C*n) 


2.5-8 


2.5-9 


The "chordwise” portion, g (£,g) , is composed of a pressure term 
from each side edge, 

g(€*n) = g (c»ri*y 0 ) -gUrneYi) + s f [ g(£rn»-y 0 ) -gUfrw-yi) 3 2.5-10 


where 


y ,y. = Outboard and inboard side edges for right hand side of 
° 1 planform 

— y , — y , = Outboard and inboard side edges for left hand side of 
1 planform 

S- = +1 for a symmetric analysis, -1 for an antisymmetric 
f analysis 


The terms g(£,n,y ) consist of a portion derived, reference 1, 
in an asymptotic ixpansion process to satisfy the change in 
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boundary conditions across the hinge and side edge, and 
modification functions which maintain the necessary singular 
characteristics at the hinge and side edge hut cause the total 
expression to have the correct characteristics at the planform 
boundaries. Two boundary value problems are used: (1) partial 

chord control, (2) full chord control. The partial chord 
expression is used for all side edges associated with trailing 
edge controls. The partial chord expression is subtracted from 
the full chord expression for all side edges associated with 
leading edge controls. 

The spanwise modification function used for both the partial 
chord and full chord expression is 

H(h) '■ * + - 5C + .375C2) , C = ln-y |/|±s-y | 2.5-11 

s 

The following coefficients are used independently of side edge 
for the partial chord control expression. 

Ci *■ +' M*/e| + -5 )/B h - i(C-C c ) (2k + 


C 2 = [ (C-5 c )^ 2 (1 -5tanA H ) - 2ik] 


2.5-13 


The contribution to g (£,n) for each partial chord side edge is 
then. 


g(S#ruy ) = H(n)[c,E t M( 5 -x ,. n -yj + c 2 E 2 < n - y jn(£-x , -y ) ] 2 . 5-14 

w 5 • b b b 


C 

t 

c 


2.5-15 


where U-^) / ( C c "C 1 ) for 

S c H t )/( 5 t - 5 c ) for 5 > E( _ 
e 2 = [ (5— C x ) 2 (5 t — S> 2 /( (C-g) ! +(n-y s ) 2 H (C t -C) 2 + (n-y s ) ! )1 1/4 


M(5-x s ,n-y s ) = 1n[ R- (B 2 (rry s ) + (?-x g ) tanA H ) /B H ] 

N(£-*V n " y s> = 1n f R "f5' x s )l 


2 . 5 - 16 

2.5- 17 

2.5- 18 


and r = /(£-x ) 2 + B 2 ' (n~y J 2 

S 2 


2.5-19 


Note that C 1 ,C 2 ,M, and N were derived using the asymptotic 
expansion process; E t and E 2 are chordwise modification functions. 

The following coefficients which are independent of (£,n) are 
side edge dependent in the full chord control expression. 
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where C t = 1 - ik(x^-x^) 


C 2 = k (- (x c ~x -j) kM 2 / B 2 + i ( 6 2 -M 2 ) /g 2 ) 2.5 

(x c - Xl ) = u c (y s ) - c 1 (y s )) 

The following terms are dependent upon side edge and upon ( 5, 0 ) 
in the full chord control expression. 


2.5- 20 

2.5- 21 


G IS = E 3 Qi/7Fq 

g^-e.^q, 

G s2 » (CiL 2 + C 2 Li) 


where 


Qi = Rsign(n~y )Li + tanA L 2 

S Li 


l 2 = 


L = 


R = 


’C J S 


’C - 1 S 


2.5-22 


2.5-23 


2.5-24 


2.5- 25 

2.5- 26 


2 . 5-27 


’c J s 


2.5-28 


and Cj = tanA T + ik[ 2 (?- ? c (y g ) ) tanA T - (3 Bf-2 B 2 ) (rryj l^B 2 


C 2 = 1 ♦ ik 2 )/2B 2 

and the chordwise modification function E 3 is. 


2.5- 29 

2.5- 30 


E 3 = /3 - 2£ - £ 2 /2 
additionally, 

G L 1 = (n-y s ) [ e 3 In ( (Cl +c 2 ) 2 ) - Evince 2 (n -y s ) 2 ,] 
G l2 = [1 + ik 2 (^-^ 1 (y s ) - .75tanA L (n-y s )) /B 2 3 g l1 
where C x = L 2 - tanA L (g-y g ) 

C 2 = n /FT7T7T L, 


2.5-33 


2.5-35 


~ ^ Z 

and the chordwise modification function E 2 is 

e = [ (C t -5) 2 { (C t -^) 2 +B 2 (n-y s ) 2 }/U t -£i> H U t ~V 2+ B 2 (n-y s ) 2 }] 1/4 

2 . 5-36 

and finally. 


G AT = E 2 Bsign (n-y g ) arctan (C 4 /C 2 ) 
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where C t = <* s > L, 2.5-38 

C 2 - Lf + (n-y g ) tanA L 2.5-39 

Ncrte that all terms except E 2 and E 3 were derived using the 
asymptotic expansion process. 

The contribution of each full chord side edge to g(5, n) is 

g( ^^y S ) = CC IS G IS +C S1 G S1 +C S2 G S2 +C L1 G L1 +C L2 G L2 +C AT G AT 3 2 -5-40 


• [e ik2fl2 { (y s ) VB 2 ] 




2.6 DOWNWASH DEFINITION 

The left hand side of the downwash integral equation is the 
kinematic downwash, or kinematic angle of attack, W (x,y) . The 
kinematic downwash, which is derived from the modal displacements 
for some structural vibration mode, is the boundary condition 
which must be satisfied by the as yet unknown pressure 
distribution under the integral equation 2.3-1. 

The following is applicable for any mode j; the subscript is 
omitted.. Reference 1 presents a more detailed derivation and 
explanation. 

If the equation of the surface of a general body in a flow field 
is written, F(x,y,z,t) = 0, the condition of no flow through 
the body is 

DF/Dt = 3F/3t + (3F/3x)u + (3F/3y)v + (3F/3z)w 2.6-1 

where DF/Dt is the substantial derivative with respect to time. 
When the body is a flat plate undergoing sinusoidal motion, 

z (x,y,t) = Z (x,y) e iwt 2.6-2 

The velocity normal to the surface, w, can be written, 

w = -[ 3z/3t + (3z/3x) u + (3z/9y) v3e la>t 

= -[(3Z/3x)u + (3Z/3y)v + iwZ]e lu)t 2.6-3 

Assuming there is no spanwise .f lew, the kinematic downwash 
(amplitude ratio) , w = (w/V) e lwt , is 

W(x,y) = -[ (3Z (x,y)/3x) (u/V) + ikZ (x,y) ] 2.6-4 

where the term (u/V), called the Velocity Pr ofile , is identically 
one for a flat plate. A first order approximation of thickness 
effects may be introduced using a velocity profile which is not 
identically one, see reference 2. This modification of the real 
part of the boundary condition is particularly significant when 
attempting to calculate control hinge moments for non-flat plate 
airfoil sections. 

The RHOIV program will, at the user's option, generate an 
additional kinematic downwash column for a gust analysis. The 
forms available are 

W (x,y) = [cos ($) ~ i»sin($) ] 2.6-5 

$ - [k(x-x^ e ^)], x re £ = a zero phase gust reference point 
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which is referred to as a gradual penetration gust, and 
W (x , y ) = [ 1 - oi] 

which is referred to as a discrete gust. 


2.7 SOLUTION FOR UNDETERMINED COEFFICIENTS 

In order to solve for the unknown coefficients of the assumed 
main surface pressure terms, a (m, j) , the kinematic downwash, 

W(x,y; j) , is first modified by removing the mathematical dcwnwash, 
W*(x,y,n;j), associated with each_control surface in the analysis. 
The resulting residual downwash, W(x,y;j) , which is smooth and 
continuous, is used in the set of linear equations which is 
solved for the unknown coefficients'. 

[W(x,y;j) ] = [W(x,y; j) ] - [W* (x,y,n; j) ] 2.7-1 

where 

[W*(x,y,n; j) ] = [C* (x,y,n;m) ][A(n;m, j) ] 2.7-2 

then [C (x,y ;m) ][a{m, j) ] = [W(x,y; j) ] 2.7-3 


2.8 UNSTEADY PRESSURES, SECTIONAL AND TOTAL GENERALIZED FORCES 

The final results generated by the RHOIV program, consist of delta 
pressures and generalized forces. The pressures are determined 
by evaluating the various pressure terms used at the desired 
output points and combining with the required coefficients. 

The generalized forces are determined by integrating pressures 
times modal displacements for all combinations of modes used. 
Sectional forces involve integrating along a chord; total forces 
involve integrating over the area of the planform half span. ' 

All program output has the coefficient of q = 0.5PV 2 . 


The pressure at any point (£,h) for some mode j is composed of 
a contribution from the assumed main surface pressure terms and 
a contribution from the pressure terms associated with each 
control surface. 


AP(S,n;j) = ap (5,n; j) + eap (?,n 

n 

where 


#n; j). 


AP r (Crn;j) = 4pV2/s2- n 2{Ap{C,n;m)} • {a (m, j) } 

= contribution from assumed main surface 
pressure terms 


2 . 8-1 

2 . 8-2 


AP c (C,n,n; j) = 4 pV 2 /S 2 -n 2 { Ap ( l, n, n ; m) } • {A(n;m,j)} 2.8-3 

S 

= contribution from control surface n 
pressure terms • 


Note that the terms Ap(£, n;m) can be calculated independently 
of k value and Mach number. The program output for pressures, 

AP (f,n; j )/q, has dimensions of (modal displacement units) /planform 
length units) . 


The sectional forces at a spanwise station n for the combination 
of i displacement mode and j pressure mode are also composed 
of contributions from main surface and control surface pressure 
terms. 

s S t (n) 

Q (n;i#j) = f Z(C,n;i) AP(?,n; j)d? 2.8-4 

Ci(n) 

[Q S (n;i,j) ] = CC®(n;i,j)] + l[ Q® (n,n; i, j) ] 2.8-5 


where 


s 


[Q r (h;i,j) ] = [ Q { 0 ; i/ m) ][ a (m, j) ] 


2 . 8-6 


= contribution from assumed train surface 
pressure terms 


[Q s (n;i,m) = 4 pV 2 /s 2 - n 2 fy ( £, n ; i) Ap ( £, n;,m) <H 

(n) 


2.8-7 


Note that Q s may be calculated independently of k value and Mach No 


s , C t (rv> _ 

(n,n;i, j) ] = 4 pV 2 /S 2- n 2 [ / Z(£,n;i) AP-(5rn»n;m) d£][A (n;m, j) ] 2.8-8 

. 3 q(n) 

= contribution from control surface n 
pressure terms 

The program output for sectional forces, [ Q s (n;ir j) ]/<3» has 
dimensions of (modal displacement units) 2 . 


Similarly, the total generalized forces are given by 

s c t ( n ) 

Q (i# j) = / / z (C,n;i) AP(c,n: j)dcdn 2.8-9 

0 h (r]) 

C:Q(i»j)] = CQ_(i»j)'} + E[Q (n;i, j) ] 2.8-10 

. : n s 

where 

[ Q r (i# jj ] = [ Q r (i»ir) ][ a (m, j) ] 2.8-11 

= contribution from assumed main surface 
pressure terms 

s C. (n) 

U (i,m) = 4 P V 2 / /V Z (5, n ;i) Ap(C/n;ni)dCdn 2.8-12 

0 q(n) 

Note that ■ may be calculated independently of k value and Mach No. 


S (n) 

[Q (n;i,j)]= 4 pV 2 [/ / x VS 2 -n i Z (£,n;i) Ap (£ ,n ,n;m) d£dr|] [A(n;m, j) ] 2 . 8-13 
S 0 ? 1 (n) 

= contribution from control surface n 
pressure terms. 
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3.0 COMPUTER PROGRAM USAGE 

3.1 MACHINE REQUIREMENTS 

The RHO IV program system is written for the CDC 6000 series 
computer. It requires the use of a card reader, line printer, 
disk storage, a minimum of zero and a maximum of five tape drives. 

3. 2 OPERATING SYSTEM ‘ ' 

The program runs under the SCOPE 3.1 cr KRONOS operating systems. 
All system routines used are assumed to be standard CCC release. 
With the exception of COMPASS routines used for shifting and 
vector inner products, all source routines are coded in CDC 6600 
FORTRAN IV. The overlay loading feature is used. 

3.3 STORAGE ALLOCATION 

The RHO IV program will load under a field length of 37000 a . 

The program has been written to use blank common as a working 
area for those portions of the analysis which are dependent upon 
the sire of the user's problem. Specifically , there is no program 
limitation on the number of modes, modal input points, pressure 
output points, or sectional force output chords. 

The minimum core requirement of 53000 a words is determined by 
C-roatrix calculation in which variable dimensioning is not used. 

The requirements for the other sections will he calculated and 
printed at execution time; they may be precalculated using the 
formulae below. Some examples of required field lengths for 
selected problems is given following the equations. Note only 
those sections which will be used need be considered. 

The following variables are used in describing core requirements; 

1. NOVP Number of user supplied velocity profiles 

2. NVPP (I) Number of points associated with velocity profile I. 

(Refer to user input, page 64) 

3. NPRC Number rressure output chords 

4. NPPT NPRC*NP?RC 

5- NSGFC Number of sectional force output chords 

(Refer to user input, page 66. Note that default 
values may be supplied by the program.) 

6. NOCS Number of control surfaces 

NIPTS(I) Number of modal input points for zone I 
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vO 00 


11 . 

12 . 

13. 

14 . 

15. 
16 

17. 

18. 
19. 

(D 



( 2 ) 


(3) 


m 


Max [NIPTS (I) , 1=1 , NOCS + 1 ] ; ; ; 

dumber of modal input, zones , =NOCS+1 
Number of displacements modes •; 

Number of pressure modes 
= NDMDS if .NOT. GUST 
= NDMDS +1 if GUST 
(Refer to user input, page 68) 

Number of downwash points = NDNC*NPD'WC 
Number of assumed main surface pressure terms 
(=NDWP currently) 

(Refer to user input, page 66) 

Length of IIA information 
Length of velocity profile information 
Length of control rotation information 
Maximum number of points/chcrd required 
for sectional or total generalized forces 
Maximum number of chords required for generalized 
force integration 
Length of CMFILE index 
(Refer to discussion below) 

The interpolation information arrays are used to calculate 
control rotation coefficients, basic downwash matrix, and 
sectional and total generalized force precalculated 
information. The length involved is, 

(a) LIIA = The sum of the lengths of all IIA arrays input 

by user (IIAIN=.TRUE.) , or 

(b) LIIA = NZONES* (3*NDMDS + 23) + (NDMDS+4) *Sum[ NIPTS (I) , 

1=1 , NZONES ] 

If controls are present, control rotation coefficients are 
used in all preparation routines. The length involved is, 

LCCR = 4*NOCS*NDMDS 

If velocity profiles are used, cubic splines are generated 
for each profile, and the information used in calculation 
of control rotations and the basic downwash matrix. 

LVP = 5*Sum[NVPP (I) ,1=1, NO VP] - NOVP 

For the purpose of performing sectional and total generalized 
force integration, maximums can be placed on the number 
of integration chords and points/chord which will be 
required. The maximum number of points/chcrd required for 
sectional or total generalized force integration is. 


MIPTS 

NZONES 

NDMDS 

NPMDS 


NDWP 

NPTRM 


LIIA 

LVP 

LCCR 

MPCHD 

MICHD 

. LINDEX 



MPCHC = 22 + (4*NOCS+12) if KGCS>0 


The maximum number of chords required for total generalized 
force integration is, 

MICHD = 17 + 8*N0CS 

(5) If a user supplied CMFILE is present, a CMFILE index is 
required. The length of a CMElLE index is, 

LINDEX = 1 3*No. Main Surface C-matrices + No. Control 
Surface C-matrices 

A minimum for each of the sections is given, followed by additive 
amounts for each subsection. The cere required is the maximum 
required for any section used. The requirement for a section 
is the section minimum plus maximum required for any subsection 
used. 


I. 


INPUT PREP APATI 0 N 

A. 


B. 


Velocity profile 
input 

Pressure results 
input 


C. Sectional force 
results input 

D. Modal input 


51000 a + Minimum of 1440 a 

7 *Max[ NVPP (I ) ,I=1,N0VP]-4 

NPRC + NPPT 

NSGFC 


la.MIPTS* (NDMDS+4) 
or b. Maximum size of user input IIA 


II. MODAL INPUT PREPA RATIO N 37000 a + LVP 

A. Calculation of inter- 1. MIPTS* (NDMDS + 2) +9*NDMDS 
polation information +2a . (MIPTS+3 ) **2MIPTS+3 

or b.23+3* (MIPTS+NDMDS) +MIPTS* 
NDMDS 


B. Calculation of 
control rotation 
coefficients 

III. RESUL T PREPARATI ON 

A. Formation of basic 
kinematic downwash 
matrix 


NZONES+LIIA+LCCR+3* NDMDS 

370 00 „ +NZONES +LIIA + LCCR 
LVP+NDWP* (2*NPMDS+3) +NOCS+1 
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B. Unsteady pressure 1. NPRC+NPPT 

results preparation +2a .NPPT*NPTRM 

or b . NPRC*NPMDS*NOCS 

C. Sectional force 1. NZ0NES+3*MPCHD+NPMDS 

preparation +2a . NDMDS*NPTRM 

or b.MFCHD*NPMDS 
+ 3. NSGFC 

D. Generalized force 1. ,2 . as in C above 

results preparation +3. 2*MICHD 

IV. C-MAT RIX CALCULATION 53000 a 

V. C-MAT RIX LIERARY USAGE 35000 a + LINDEX + NDWP*NPTRM*2 

VI. SOLUTION SECT ION 35000 a 

A. C-matrix/downwash la. 2*NDWP*NPTRM 

matrix printed or b.NOCS* (8*NDWP+4*NDMDS) +8*NDWP 

output 

B. Solution for 1. NDWP* (2*NPTRM+N) 

coefficients of main +2 . ,NOCS* (8*NDWP+4*NPMDS) 
surface pressure 

terms 

Note that the value of N is selected on the user's field 
length, N=minimum of 3, maximum of 1+2*NPMDS. 

C. Pressure coefficient 8*NPTRM 
output 

VII. RESULTS 35000 a 

A. Unsteady pressure 1. NPRC+NPPT* (NPTRM+6) 

calculation +2a. 2*NPTRM 

or b. 2*NFRC*N0CS 

or C. 4 *NFFT 

+ 3. 2*NPPT (if NOCS>0) 

B. Sectional force 1. 2*NDMDS*NPMDS 

calculation +2a. NDMDS*NPTRM+2*NPTRM 

or b. 119 + 20*NOCS (if NOCS>0) 

+3. NSGFC 

C. Generalized force 1 . , 2 - as in B above 

calculation 



In order to provide a user with an easily determined initial 
field length estimate the preceding equations have been applied 
for several combinations of user controlled parameters. 

The following field length requirements are for a lifting surface 
with two (2) control surfaces where the assumption is made that 
the maximum number of input points per input zone is no more 
than 75% of the total number of input points. The field lengths 
have been rounded uo to the nearest 5000 a . 


Total No. 

No. 

No. Eownwash 

No. Pressure 

FL Required 

Input Pts . 

Modes 

Pts. 

Output Pts. 

1000 a 


72 

231 


100-0 

35 

231-100-0 

72-35 

231-100-0 

72-35 

231-100-0 



150 


150 


130 


115 


125 


105 


110 


65 


70 


60 


105 


60 


65 


minimum 



231-100-0 


231-100-0 


231 


100-0 


231 


100-0 


231 


100-0 


231 


100-0 


231 

105 

100-0 

60 

231 

65 

100-0 

minimum 



























































3.4 TIMING 

The central processor time, required to execute any problem is 
almost entirely dependent upon the number of C- matrices which 
must be calculated. A large (>90%) reduction in CP (central 
processor) time may be gained by utilizing previously generated 
C -matrices (refer to Section 3.5.2). 

The principle factors in the CP time required to calculate a 
C-matrix are (1) the number of downwash points, (2) whether the 
matrix is associated with the main surface, a trailing edge 
control surface, or leading edge control surface, and (3) whether 
the condition is steady state (k=0) , Or unsteady (k*0) . 

Items (2) and (3) plus planform shape and position of a downwash 
point on the planform will determine the CP time required for 
each downwash point. The values given below are average 
requirements for the sample problem of Section 3.9. 

SURFACE CO NDITION CP SECONDS/DOWN WASH PT .' 


MAIN 

k 

k 

II A 

0 

0 

1.0 

3.0 

TRAILING EDGE 

k 

— 

0 

3.0 

CONTROL 

k 

> 

o 

11.0 

LEADING EDGE 

k 

zz 

0 

11.0 

CONTROL 

k 

> 

0 

27.0 


The above values may be used initially to estimate CP time 
required fcr a users problem. Note that a C-matrix must be 
calculated or retrieved from a library for the main surface and 
each control surface in the problem for each condition. At 
execution time the RHO IV program prints the specific CP time 
required per downwash point for any C-matrix calculated as well 
as printing a breakdown of CP usage by the various other sections 
of the program. 



3.5 FILE I/O 

3.5.1 Fil e Utili za tio n 

RHO IV uses standard input and output, two internal scratch 
files, and up to five (5) user specified input/output files. 

The scratch files are referenced in the program by different 
names, dependent upon the usage. All user input/output files 
have the record format described below for REALM/ WRIT EM. 

The files referenced are: 

(1) INPUT Standard input (BCD) 

(2) OUTPUT Standard output (BCD) 

(3) RH0SC1 Scratch file (binary) 

(MISFILE) Modal input scratch file 
(RESFILE) Result scratch file 

(4) RH0SC2 Scratch file (binary) 

(INSFILE) Input scratch file 
(CMS FILE) C-matrix scratch file 
(COFFILE) Coefficient file 

(5) MIFILE Modal input file, user specified file for 

input modal displacements and associated 
points (binary) 

(6) CMFILE C-matrix file, user specified input/output 

file containing the library of previously 
calculated C- matrices, (binary) 

(7) DPFILE Pressure output file, user specified 

output file containing all unsteady 
pressure results 

(8) SGFFILE Sectional force output file, user specified 

output file containing all sectional force 
results . 

(9) GFFILE Generalized force output file, user specified 

output file containing all generalized force 
results . 

Note that the user specifies the file names for (5) - (9) in data 
input. Files (5) , (7) - (9) may be equivalenced in any combination 
using one or more file names, the output for (7) -(9) will then 
be inter leafed on a k -value, Mach number condition basis. 


If a CMFILE is defined it must be descrete. 

The user may also specify an initial file position for (5) - (9) . 
MIFILE is positioned and used, then DPFILE, SGFFILE, and GFFILE 
are positioned, in tha t order . Therefore if a user gives two 
or more files the same name, the initial file position specified 
for the last in the above sequence would be used. 

All input/output files use a two record format for each array 
written. Routines READM and WRITEM are used to read and write 
these Records . 

READM/ WRIT EM FOR MAT 

Record 1 - ID Record Word 1 - 5HMATI0 

2 - MR0W = Row length of array 

3 - MC0L = Column length of array 

4 - LID = Length of user ID 
5- (4+LID) - ID = User ID 

Record 2 - ARRAY Record Words MR0W*MC0L - ARRAY, written 

(ARRAY (I, J) ,I=1,MR0W) ,J=1,MC0L) 

3.5.2 C-Matrix Librar y 

The C-matrix library, maintained on a user specified 1/0 file, 
CMFILE, consists of a two file set. File 1 consists of all saved 
C-matrices written sequentially in the order saved- File 2 
consists of an index of File 1. If the user specifies a CMFILE, 
the index is examined for a match each time a C-matrix is 
required. If a match is found the desired C-matrix is accessed 
and used. If a match is not found, the desired C-matrix is 
calculated, the index updated, and the new C-raatrix and updated 
index are written to CMFILE. If a legitimate index is not found 
on the first attempt to access information in an execution, the 
specified file is assumed to be a new file; an initial null index 
will be generated and execution will continue. There is no 
program limit on the number of C-matrices in the library. 

All arrays on CMFILE are written and read using READM/WRITEM. 

The CMFILE index consists of a list of entries associated with 
each C-matrix sufficient to allow testing for a match on a 
retrieval attempt. Note that the user must utilize MSID in input 
to deferentiate between main surface planforms, and CSID to 
deferentiate between control surfaces associated with a main 
surface planform. A main surface entry in the INDEX consists 
of 13 words, and contains a counter for the number of associated 
control surface entries which follow. A control surface entry 


consists of one word; it can only be referenced through its 
associated main surface entry. 

MAIN SURFACE ENTRY , • 


WORD 

EITS 

.VAR., TYPE 

VARIABLE. 

DESCRIPTION • 

1 

59-18 

H 

MSID 

Main surface ID 


14-00 

I 

MATP0S 

Position of matrix in file 
CMF1 

2 

35-30 

I 

SYM 

Symmetry indicator 


29-24 

I 

NDWC 

No. downwash chords 


23-18 

I 

NPDWC 

No. points/downwash chord 


17-12 

I 

NSPT 

No. s panwise main surface 
pressure terms 


11-6 

I 

NCPT 

No. chordwise main surface 
pressure terms 


5-0 

I 

NCSE 

Number of associated con- 
trol surface entries 

3 


R 

S 

Semi-span 

4 


R 

KVAL 

uj/V - reduced frequency 

5 


R 

MACH 

M - roach number 

6 


H *■ 

DATE 

Entry date 

7-13 


H 

RTITLE 

- En t ry run t it le 


CONTROL SURFACE ENTRY 


WORD BITS VAR. TYPE I VARIABLE 


DESCRIPTION 


59-18 

17-15 


14-00 


CSID 

CSTYPE 


MATP0S 


Control surface ID 
Control surface type 
Full span y 
Tip span I Leading edge 

Mid span I or 

Partial span' Trailing edge 
Position of matrix in 
file CMF1 


New mam surface entries are appended to the bottom of the 
index. New control surface entries are inserted fcplow the last 
previous control surface entry of its associated main surface 




















3.5.3 File formats 


All files, with the exception of INPUT and OUTPUT are described 
in this section; the order is alphabetical. 

A file name, equivalent name, file type, and short description 
of usage are given in addition to record formats. 

With respect to the headings, the following applies. 


(1) REC. N0. 

(2) REPETITION 


RECORD N0. OR IDENTIFIER 

AN INTEGER SPECIFYING THE NUMBER OF IMMEDIATE 
REPETITIONS OF THE RECORD. REPETITIONS OF SETS 
OF RECORDS ARE DESCRIBEE PRIOR TO THE SET. 


(3) 


(*0 


VARIABLE 

(DIMENSIONS) 


PROGRAM, INPUT, OR OUTPUT VARIABLE NAME. 

NUMBER OF ELEMENTS ASSOCIATED WITH THE \ 
VARIABLE. NOTE THAT ALL ARRAYS ARE WRITTEN 
IN CDC6000 STORAGE ORDER, E.G., 


A(I,J,K) 


( ( (A (i , j, k) ,i=1,I) ,j=1,J) ,k=i,K) 


TYPE OF INFORMATION 

I - INTEGER . 

R - REAL 

C - COMPLEX (NO. TERMS=2*FIRST DIMENSION) 
H - HOLLERITH 
M - MIXED 


(5) DESCRIPTION I DESCRIPTION OF INFORMATION EY VARIAELE 
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The following variable names are used in describing DIMENSIONS. 

I. N0CS NUMBER OF CONTROL SURFACES 

N0VP NUMBER OF VELOCITY PROFILES 

NDWP NUMBER OF DOWNViASH POINTS 

NPTRM ■ NUMBER OF ASSUMED MAIN SURFACE PRESSURE TERMS 
(=NSPT*NCPT, in most cases = NDWP) 

5. NDMDS NUMBER OF DISPLACEMENT MODES 

6. NPMDS NUMBER OF PRESSURE MODES (EITHER NDMDS OR NDMCS+1) 

7. N0KVAL NUMBER OF REDUCED FREQUENCIES 

8. N0MACH NUMBER OF MACH NUMBERS 

9. N0C0ND NUMBER OF CONDITIONS = N0KVAL*N0MACH 

10. NPRC NUMBER OF PRESSURE REPORT CHORES 

II. NPPRC NUMBER OF POINTS/PRESSURE REPORT CHORD 

12. NPPT NUMBER OF PRESSURE OUTPUT POINTS = NPRC*NPPRC 

13. NSGFC NUMBER OF SECTIONAL FORCE OUTPUT CHORDS 

14. NPZ0NE NUMBER OF INPUT POINTS/MODAL INPUT ZONE 

15. LINDEX LENGTH OF CM FILE INDEX 










FILE NAME: CMFILE FILE TY-EEL; SEQUENTIAL BINARY 

CMFILE is a user specified INPUT/OUTPUT file generated and used 
by the RHOIV program as a library of previously c coated 
matrices. It is a multifile file; the first file is the set>- 
of all saved C-matrices, the second is. an index of the firsts. 

Routines READM/WRITEM are used to READ/WRITE the ' V- ; V-V 

REC. 

NO. 

REPETITION 

: VARIABLE (DIMENSIONS). 

T 

DESCRIPTION 

The record pair 1A,B is repeated NOCM times (refer to INDEX deec. ) .. 

1A 

i 

READM/WRITEM ID INFO (4) 
USER ID (4) 

1- 8HC-MATRIX 

2- S . - . . : ■ 

3- KVAL 

4 - MACH 

•t 

E 

R 

R 

R 

•4 

REFER TO PAGE 38 

ARRAY NAME 
SEMI- SPAN 

REDUCED FREQUENCY=W/V 
MACH NUMBER 

IB 

i 

t • . 

■ vi 

C (N,NDWP) 

C 

COMPLEX C- MATRIX, 
.SAVED IN TRANSPOSE 
FORM. N=NPTRM 
FOR A MAIN SURFACE 
C-MATRIX, AND =4 FOR 
A CONTROL SURFACE 
C-MATRIX 

2 

i 

END OF FILE 

3A 

i 

READM/WRITEM ID INFO (4) 
USER ID (5) 

1- 10HCMFL INDEX 

2 - NOCM 

3- NMSNTRY 

4- CDATE 

5- LMDATE 

1 

, 

, ■ 

H 

I 

I 

H 

H 

REFER TO PAGE 38 

ARRAY NAME 
NUMBER OF C-MATRICES 
NUMEER OF MAIN SURFACE 
ENTRIES OR C-MATRICES 
LIBRARY CREATION DATE 
LAST MODIFICATION DATE 

3B 

i 

INDEX (LINDEX) 

M 

CMFILE INDEX - REFER 
TO PAGE 38 
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FILE NAME: CMSFILE {= RH0SC2) FILE TYPE: SEQUENTIAL BINARY 

CMS FILE is a sequential binary file used internally lay RHOIC. 

It is generated by CMCALC and/or RDWRTC and used by PCMDWM and ... 
SCLUTON during the k- value, Mach no. condition cycle. ‘ 

REC. 

NO;” 1 ' 


VARIABLE ( E I ME N S I O.N£) 

"X 

T 

DESCRIPTION 

■ 

NDWP 

( ROWC(NPTRM) 

1 

ROW OF MAIN SURFACE 
PRESSURE TERM C -MATRIX 

IF NOCS>0, a set of Rec. 2 is repeated per control, i.e. NOCS times. 

2 

ndwp 

RONC (4) 

B 

ROW OF CONTROL SURFACE 
PRESSURE TERM C-MATRIX 
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FILE NAME: 


COFFILE (= R HO SC 2) 


FILE TYPE: SEQUENTIAL BINARY 


COFFILE is a sequential binary file used internally by 'RHO.IV. 
It is generated by SOLUTON and used by PCMSPT and the various 
result routines, FORMDP, FORMQSy and FORMQ , during the k -value 
Mach no. condition cycle. 


REC. REPETITION 
MG. 


VARIABLE (DIMENSIONS) 


DESCRIPTION 


NPMDS 


CMSPT (NPTRM) 


COEFFICIENTS OF MAIN 
SURFACE. PRESSURE TERMS 
FOR A PRESSURE MODE 











FILE NAME: DPFILE , FILE TYPE: SEQUENTIAL BINARY 


DPFILE is a user specified output file which will contain any 
unsteady pressure results. Routine WRITEM is used to produce 
an ID record and record of values. 


EBB 

REPETITION 

VARIAELE (DIMENSIONS) 

a 

DESCRIPTION 

1A 

1 

R EADM/WRITEM ID INFO ( 4) 
USER ID (4) 

1- lOHRHOIV Y-DP 

2- NPRC 

3- NPPRC 

4 - NPPT 

H 

I 

I 

I 

REFER TO PAGE 38 

ARRAY NAME 

NUMEER OF PRESSURE 

REPORT CHORDS 

NO. POINTS/PRESSURE ; 

REPORT CHORD 

NO. PRESSURE POINTS ' 

IB 

1 

YPRC(NPRC) 

R 

PRESSURE. REPORT CHORDS 

2A 

1 

READM/WRITEM ID INFO (4) 
USER ID (4) 

1- 10HRHOTV X-DP 

2- NPRC 

3- NPPRC 

4 - NPPT 

I 

H 

I 

I 

I 

REFER TO PAGE 38 
AS ABOVE 

2B 

1 

XPPT(NPPT) 

R 

PRESSURE REPORT POINTS 


The record pair 3A,B is repeated per mode, i.e. NPMDS times, 
for each k-value, Mach no. condition. The record pair will occur 


a total of NO KVAL*NOMACH* NPMDS times. 


3A 


1 


READM/WRITEM IE INFO 
USER ID (6) 

1- 8HRH0IV DP 

2- RKVAL 

3- BO 


4-S 


5 -MACH 


6-IMD 



REFER TO PAGE 


38 


H 

R 

R 

R 

R 

I 


ARRAY NAME 

REF. K~VALUE= bow/V 

K-VALUE REF. LENGTH 

SEMI-SPAN 

MACH NUMBER . 

MODE NUMEER 


3B 


1 


PRESS (NPPT) 


COMPLEX PRESSURE AT 
OUTPUT POINTS FOR MODE 
IMD. 


C 






























FILE NAME: GFFILE , , FILE TYPE: SEQUENTIAL BINARY 

GFFILE is a user specified output file which will contain any 
generalized force results. Routine WRITEM is used to produce 
an ID record and a record of values. 


REC. 

NO. 


REPETITION 


VARIABLE (DIMENSIONS) 


DESCRIPTION 


Record pair 1A,B are repeated per k-value, Mach number condition. 
The record pair will occur a total of NOKVAL*NOMAOH times. 


1A 

1 

READM/WRITEM ID INFO 

(4) 

M 

REFER TO PAGE 42 



USER ID (5) 
1-7HRH0IV Q 


H 

ARRAY NAME 



2-RKVAL 


R 

REFERENCE K-VALUE = 
bo oj/v 



3-BO 


R 

K-VALUE REF. LENGTH . 



4-s 


R 

SEMI-SPAN 



5 -MACH 


R 

MACH NUMBER 

IB 

1 

Q (NDMDS, NPMDS) 

C 

COMPLEX GENERALIZED'- 
FORCE MATRIX. 


















FILE 

NAME: INSFILE (= RHOSC2) FILE 

TYPE:' SEQUENTIAL BINARY 

INSFILE is a sequential binary file used internally by RHOIV. 

It is generated by various input routines and used by the various 

preparation routines IIACAL, CRCOEFF, RESPREP, 

VPRDR . 1 

REC. 

NO. 

r— : — 

REPETITION 

VARIABLE (DIMENSIONS) 

T 

/DESCRIPTION 

Record set 1 occurs if NOVP>0, and is repeated NOVP times. 

1A 

1 

YVP 

R 

STATION AT WHICH A 
VELOCITY PROFILE IS 
DEFINED 



NVPP 

I 

NUMBER OF POINTS AT 
WHICH VELOCITY PROFILE 
IS SPECIFIED. 

IB 

1 

XVP (NVPP) 

R 

VELOCITY PROFILE POINTS 



CVP (4 ,NVPP-1) 

R 

CUBIC COEFFICIENTS IN 
XVP DEFINING PROFILE 

Record set 2 occurs if any pressure results 

have been requested. 

i.e. 

if DPPRT^O, or DPFILE^O. 



2A 

1 

YPRC(NPRC) • 

R 

PRESSURE REPORT CHORDS 

2B 

1 

XPPT(NPPT) 

R 

PRESSURE REPORT POINTS 

Record set 3 occurs if any sectional force results have been 

requested, i.e. if 

SGFPRT^O, or SGFFILE^O. 



3 

1 

YSGFC (NSGFC) 

R 

SECTIONAL FORCE OUTPUT 
CHORDS 

Record set 4 occurs unless IIAIN=.TRUE. , i.e 

. unless interpolation 

information arrays 

have been input directly. 

The record set 

is repeated once for each input zone, i.e. NOCS+1 times. 

4 A 

1 

NPZONE 

I 

NUMBER OF INPUT POINTS 
IN ZONE 



X (NPZONE) 

R 

X VALUES OF INPUT 
POINTS 



Y (NPZONE) 

R 

Y VALUES OF INPUT 
POINTS 

4B 

NDMDS 

Z (NPZONE) 

R 

MODAL DISPLACEMENTS AT 
INPUT POINTS FOR ONE 
MODE 
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FILE NAME: MIFILE FILE TYPE: SEQUENTIAL BINARY 

MIFILE is a user specified input file which is used to provide 
modal input points and displacememts to R HO IV from disk, or tape. 
Routine WRITEM, or an equivalent should be used to generate the 
ID record and value record. - 


aa 

REPETITION 

VARIABLE (DIMENSIONS) 

i 

DESCRIPTION. 

E33I 



■ 



MIFILE may consist of either records 1A,B and 2A,E or 3A, B. 

If IIAlN=.T. in input, records 3A,B should be used. otherwise 
use 1A,B,2A,B. j ‘h-, " " V V 

The set 1A,B 2A,B may occur once or 1+NOCS times. In the latter 
case, assuming NOCS/0 , the modal input points and displacements 
are associated with the main surface, and control surfaces in 
the order defined in input. 

1A 

i 

READM/WRITEM ID INFO (4) 
USER ID (10) 


REFER TO PAGE 38 . : : 

A USER ID OF UP TO 10 
WORDS. MAY BE INCLUDED. 
THE ID IS NOT USED BY 
RHOIV. 

IB 

i 

XY (NPZONE, 2) 

R 

MODAL INPUT POINTS 

2A 

i 

READM/WRITEM ID INFO (4) 
USER ID (10) 

M 

, REFER TO PAGE 38 

A USER ID OF UP TO 10 
WORDS MAY BE INCLUDED. 
THE ID IS NOT USED BY 
RHOIV. 

2B 

i 

Z (NPZONE, NDMDS) 

R 

MODAL DISPLACEMENTS AT 
INPUT POINTS 


NOTE: If there is only one input zone, i.e. NOCS=0, or if 

points for all zones are being input in the same block, 
NPZONE should be the total number input points. 
Otherwise NPZONE should be the number of input points 
for the particular zone. 


Set 3A,B is used if IIAIN=.T., in which case the pair 3A,B should 
be repeated for each input zone, i.e 1+NOCS times. 


3A 

1 

READM/WRITEM ID INFO (4) 

M 

REFER TO PAGE 38 



USER ID (10) 


A USER ID OF UP TO 10 
WORDS MAY BE INCLUDED 
THE ID IS NOT USED BY 
RHOIV. 






















FILE NAME: 


MISFILE (=RHOSCl) 


FILE TYPE: SEQUENTIAL BINARY 


MISFILE is a sequential binary file used internally by RHO IV. 
It is generated by MIINCK or IIACAL and read by IIARDR in the 
preparation prior to the k-value f Mach no. cycle. 

MISFILE will be used if any modes exist. 


REPETITION 

VARIABLE (DIMENSIONS) 

NOCS+1 

NIIA 

IIA (N IIA) 


DESCRIPTION 


NUMBER OF IIA ELEMENT 
INTERPOLATION INFOR. . 
ARRAY FOR AN INPUT ZONE 










FILE 

NAME: RESFILE (=RHOSCl) FILE 

TYPE: SEQUENTIAL BINARY 

RESFILE is a sequential binary file used internally by RHOIV. 

It is generated by 

the various preparation routines, DWPREP, 

PRSPREP, SGFPREP, and GFPREP, prior to the k 

-value, Mach no. 

cycle, and used by the solution and result routines, PCMDWM, 

SOLUTON, FORMDP, FORMQS , and FORMQ, during the 

condition cycle. 

REC. 

NO. 

REPETITION 

VARIABLE {DIMENSIONS) 

T 

DESCRIPTION 

1A 

NOCS 

IF NO CSX) 

GCR (4 , NPMDS) 

R 

CUBIC COEFFICIENTS OF 
CONTROL ROTATION 

IB 

NPMDS 

W (NDWP) 

R 

COLUMN OF THE BASIC 
DOWNWASH MATRIX 

Record set 2 occurs if any pressure results have been requested. 

i.e. 

if DPPRT/0 or 

DPFILE/0. 



2A 

1 

YPRC(NPRC) 

R 

PRESSURE REPORT CHORDS 

2B 

1 

XPPT(NPPT) 

R 

PRESSURE REPORT POINTS 

2C 

1 

XBAR(NPPT) 

R 

NON-DIMENSIONAL (BAR 
NOTATION) REPR. 

OF PRESSURE REPORT 
POINTS 

2D 

1 

MSPTRM (NPFT , NPTRM) 

R 

ASSUMED MAIN SURFACE 
PRESSURE TERMS EVAL. 
AT PRESSURE REPORT 
POINTS. 

2E 

NPMDS 
IF NOCS>0 

FETA (NPRC, NOCS) 

R 

ROTATION OF CONTROL 
SURFACES AT PRESSURE 
REPORT CHORDS FOR A 
MODE 

Record set 3 occurs if any sectional force results have been 
requested, i.e. if SGFPRT/0 or SGFFILE/0. 

3A 

1 

YSGFC (NSGFC) 

R 

SECTIONAL FORCE OUTPUT 
CHORDS 

(1) The sequence 3B-3E is repeated per sectional force output 
chord, i.e. NSGFC times. 

3B 

1 

REFER TO (1) 

QSMSPT (N DMD S , NPTRM ) 

R 

INTEGRALS OF ASSUMED 
MAIN SURFACE PRESSURE 




TERMS TIMES MODAL 
DISPLACEMENTS ALONG 
A SECTIONAL FORCE 
OUTPUT CHORD 


(2) The sequence 3C-3E is repeated per control surface if NOCS>0. 


REFER TO (1) 
AND (2) 


NIPTS 

XIPT (NIPTS) 
QTYPE (NIPTS) 


NUMBER OF QUADRATURE 
POINTS 

QUADRATURE POINTS 
QUADRATURE TYPE ASSOC, 
WITH QUADRATURE POINT 


REFER TO (1) 
AND (2) 

3E NDMDS 

REFER TO (1) 
AND (2) 


FETA (NPMDS) 


WZ (NIPTS) 


CONTROL ROTATIONS AT 
CHORD FOR ALL MODES 


QUADRATURE WEIGHTED • 
DISPLACEMENTS AT XIPT 
FOR A MODE 


Record set 4 occurs if any total generalized force results have 
been requested, i.e. if GFPRT^O, or GFFILEXO. 



QMSPT (NDMDS, NFTRM) 


NICHD 


YICHD (NICHD) 


SURFACE INTEGRALS OF 
MAIN SURFACE PRESSURE 
TERMS TIMES MODAL 
DISPLACEMENTS 

NUMEER OF SPANWISE 
QUADRATURE CHORDS FOR 
SURFACE INTEGRATION 
QUADRATURE CHORDS 


( 3 ) The sequence 4C-4E is repeated once for each control (NOCS>0) 
by quadrature chord, i.e. NOCS*NICHD times. 


REFER TO (3) 


NIPTS 

XIPT (NIPTS) 
QTYPE (NIPTS) 


NUMEER OF QUADRATURE 
POINTS 

QUADRATURE POINTS 
QUADRATURE TYPE ASSOC. 
WITH QUADRATURE POINTS 


REFER TO (3) 

4E NDMDS 

REFER TO (3) 


FETA (NPMDS) 


WZ (NIPTS) 


CONTROL ROTATIONS AT 1 
QUADRATURE CHORD 

QUADRATURE WEIGHTED 
DISPLACEMENTS AT XIPT 
FOR A MODE 




















FILE NAME: 


SGFFILE 


FILE TYPE: SEQUENTIAL BINARY 


SGFFILE is a user specified output file which will contain any 
sectional force results. Routine REAEM is used to produce an 
ID record and a record of values. 


REC. 

NO. 


REPETITION 


VARIABLE (DIMENSIONS) 


DESCRIPTION 


1A 


IB 


READM/WRITEM ID INFO 
USER ID (2) 

1 - 10HRHOIV Y-QS 

2 - NSGFC 


YSGFC (NSGFC) 


(4) 


REFER TO PAGE 38 

ARRAY NAME 

NO. SECTIONAL FORCE 

OUTPUT CHORDS 

SECTIONAL FORCE OUTPUT 
CHORDS 


The record pair 2A,B is repeated per sectional force output 
chord, i.e. NSGFC times for each k-value, Mach. no. condition. 
The record pair will occur a total of NOKVAL*NOMACH*NSGFC times. 


2A 


2B 


READM/WRITEM ID INFO (4) 
USER ID (6) 

1 - 8HRH0IV QS 

2 - RKVAL 

3 - BO 

4 - S 

5 - MACH 

6 - ICHD 

QS (NDMDS , NPMDS ) 


M 

H 

R 

R 

R 

R 

I 


REFER TO PAGE 38 

ARRAY NAME 
REFERENCE K-VALUE=bou)/V 
K-VALUE REF. LENGTH 
SEMI-SPAN 
MACH NUMBER 
OUTPUT CHORD NUMBER 


COMPLEX SECTIONAL FORCE 
MATRIX FOR CHORD ICHD.| 
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3.6 CONTROL CARDS 
There are basically four inodes of execution, from 

a. source in source form; 

b. source in UPDATE form; 

c. relocatable binary; 

d. absolute binary. 


In the following, use of specific control cards has been avoided; 
rather the required sequence of operations is specified. All 
file names with the exception of RHOIV, are arbitrary. Note 
that all overlays have the name RHOIV, thus a file RHOIV is 
generated at load time. For the cases above 

a. (1) obtain a source file, PROG (2 records, from permanent 
storage (cards, tape, permanent disk file, etc. ) . 

(2) Compile first record placing relocatable binary on 
BPROG. 

(3) Compile second record placing relocatable binary on 
BSUBS . 

(<*) Generate an alternate library on SUBLIB from BSUBS. 

(5) Load BPROG using alternate library SUBLIB. 

(6) Execute from RHOIV. 

♦NOTE; On some systems, steps 4 and 5 of a. may be combined 
into one operation, e.g., leading BPROG using BSUBS 
directly as an alternate LIBRARY. On other systems 
in which the loader has no alternate library capability, 
a preload operation may be performed in which some 
program other than the loader searches for references 
in the routines in BPROG to routines in BSUES generating 
a file LPROG, is then processed by the loader (i.e., 
the above program performs the alternate library 
function) . 


(1) Obtain an old program library file, OLDPL, from 
permanent storage. 


b 


(2) Using UPDATE generate a source file, PROG (2 records) . 
(In UPDATE terminology, PROG would correspond to the 
COMPILE file.) 

(3) Proceed with steps 2 - 6 of a. 

(1) Obtain a relocatable binary file, BPROG , main routines, 
from permanent storage. 

(2) Obtain a relocatable binary file, BSUBS, alternate 
library subroutines, from permanent storage. 

(3) Proceed with steps 4 - 6 of a. 


(1) Obtain an absolute binary file, RHOIV, from permanent 
storage. 


(2) Execute from RHOIV. 


3.7 


PROGRAM INPUT 


3.7.1 General Remarks 

The input to RHOIV consists of both BCD input, e.g. cards, and 
binary, e.g. CMFILE or MIFILE. The card input includes planform 
description, definition of user I/O files, printed output 
specifications, list of k-values and Mach numbers and modal input 
description. Card input may also include velocity profile 
definitions, the distributions of downwash points, pressure 
output points, and sectional force output chords, and modal 
displacements with associated points. The binary input may 
consist of a library of C-matrices, and modal displacements with 
associated points or interpolation information arrays. 

The card input consists of field dependent input and free field 
input. The field dependent input is identified in the field 
column of data stacking as a specific field (number) with 
associated format or as a LIST indicating sequential input per 
the FORMAT using as many cards as required. The free field input 
is identified by NAMELIST in the field column with associated 
list name in the FORMAT column. Some of the features of namelist 
input are: 

(1) Card (s) field consist of columns 2 through 80, 

(2) List consists of a $ list name in column 2 followed by a 
series of specifications continued on as many cards as 
required and terminated by a $ . 

(3) Specifications are of the form: 
a- Vname = Value 

b. Vname (1) = Valuel, Value2 , . . . , Valuen 

Where Vname is one of the variable names for the list, value 
is the associated value (s). Value may be an integer, a floating 
point number in normal or exponential form, or in the case of 
a logical variable of the form. 

. T. or .True, indicating true 

.F. or .False, indicating false 

(4) - Specifications must be separated by commas. NOTE there 

is no comma between the last specification and terminating 

$. 


(5) Embedded blanks are allowed except within the $ list name, 
variable name, or value. Note at least one blank must 
separate the $ list name and the first specification. 

(6) The order of appearance of variables on the card(s) is not 
important - the spelling is. 

(7) Any or all of the variables may be left out of the list, 
e.g., $list name .. $ is legitimate. This assumes of course 
that there is a legal default value associated with the 
variable (s) not included in the list. 

There are a number of input sets consisting of x and y locations 
on the planform. Where feasible, the option of specifying this 
information in physical or local non-dimensional coordinates 
(bar notation) has been allowed. 

3.7.2 Limitation s 

The following are size limitations within the program (also noted 
in Data Stacking) : 

2<NLE< 1 0 
2<NTE<10 
0<NOCS<6 
1 < | NDWC | <N 
1<|NPDWC|<8 
1 <NOKVAL< 2 0 
1<NOMACH<20 

The following are analysis limitations: 

No downwash chord should be placed at the tip or control surface 
side edge. In general downwash chords should satisfy |y~n |>.02S 
where n_ is the tip or a control side edge. 

o 


No. leading edge definition points 

No. trailing edge definition points 

No. control surfaces 

No. downwash chords, N=72 ./ |npdwc | 

No. Pts. per downwash chord 

No. reduced frequencies 

No. mach numbers 


It is not recommended that a downwash chord be placed at or near 
a spanwise planform discontinuity, e.g., a change in leading 
edge or trailing edge slope. 

No downwash point should be placed at the leading edge, trailing 
edge, or control surface hinge. 

The downwash point distribution should be such that the boundary 
conditions are sufficiently defined. (A downwash point on a 
control surface is not specifically required.) 

No pressure report chord or sectional generalized force report 
chord should be placed at a control surface side edge. 
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No pressure report point should be placed on a control surface 
hinge. 


DATA INPUT FLOW 


aaa 

Title 



r 


IPlanform 

Descri 

Dtion 




CARD SET 2B 


Leading Edge 
Definition 


CARD SET 2C 


Trailing Edge 
Definition 



- 






Control Surface 
Definition 





l/l 

OJ 








mmm 





< A 
<D 


CL. 

> 

o 




(Card set 3 

File 

Definii 

tion 




Printec 

Definii 

j Output 
tion 





[CARD SET 5A 
Planform Distj 
Definition 



ftARD SET 5D 


Sectional Forc4 
Output Chords 


CARD SET 6 

Condition 

Definition 


CARD SET 7A 
Modal Input 
Definition 



■ 

1 



a 


I 







— ~ 

to 

CD 

o 

o 

to 

CD 

— -f»- 

’ 

ill 


E 

•r— 

+ 

r““ . 

E 

•r— . 

[CARD SET 7C 

Pressure Output 
Pt. Definition 


o 

QC 

s- 

o 

CO 

Q 

S 

Modal 

Displacements 



§E 

CD 

a 

Q 

2: 




C/I 

CD 


</> 


CARD SET 7B 
Modal Input 
Points 


6 



FIGURE 2 





















FIGURE 2 (continued) 

; f ' _ . 
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RHOIV COSINE' DI ST RIBUTIONS 

N Distribution of chords - j = cos (in/ (2N+1H . i = 1, N 
1 .5000 


2 

.8090 

.3090 
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.9397 
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.5000 

. 1736 
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.80 90 
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.5000 

.3090 .1045 

8 
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.9325 

.8502 

.73 90 

.602 5 

.4457 .2737 .0923 

9 

.9864 

. 9458 

.8795 

.7891 

.6773 

.5469 .4017 .2455 

N 

Distribution 

of points/chord 

- Z = 

-cos f2iir/ (2N+1H , i 

1 

.5000 
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FIGURE 3 




2A 

RHOA 

NAM ELI S ' 

: NLE 

Number of leading edge defn. 
pts. 2<NLE< 1 0 (NO Default) 




NTE 

• ' ' -f 

Number of trailing edge defn. 
pts. 2<NTE<10 (No Default) 




NCOS 

Number of control surfaces 
0<NOCS<6 (Default is 0) 




NO VP 

Number of velocity profiles 
0<NOVP (Default is 0) 




MSID 

A unique identifier which will 
be associated with any main ; 
surface C -matrix generated : 

1-5 digit integer 
(Default is 0) 




SYM 

Symmetry Indicator 

0 - Symmetric (default) 

1 — Antisymmetric 




BO 

k-value reference length 
BO > 0 (Default is root semi- 
chord as determined from 
XLE (1) , XTE (1) 



■ - . ; .. 

S 

Semi-span, S> 0 (Default is 
YLE (NLE)-YLE (1) ) 

2B 

LIST 

6F10.0 

XLE (I) 
YLE (I) 

NLE pairs of X,Y values 
Defining the planform leading 
edge YLE (1+ 1 ) >YLE (I) 

2C 

LIST 

6F10.0 

XTE(I) 

YTE(I) 

NTE pairs of X,Y values 
defining the planform trailing 
edge YTE (1+1 )> (YTE(I) , 

YTE(I) = YLE (1) 






















CARD [FIELD NO 
SET COLUMNS 



DESCRIPTION 


(2) PLANFORM DEFINITION (Continued) 


If NOCS > 0, repeat set 2D NOCS times; otherwise,: omit 2D. 


1-10 A7 , 3X CSID An identifier to be associated 

With a control surface 
"C-matrix. 

11-20 F10.0 XHLI X value of inboard hinge point 

2.1-30 F10.0 YHLI - Inboard control side edge. 

31-40 FI 0.0 XHLO X value of outboard hinge 

point . ' • : 

41-50 F10.0 YHLO Outboard control side edge. 

51-60 F10.0 LEIND Leading edge control indicator 

#0 - leading edge, 

=0 - trailing edge 
(Default is 0) . 


If NOVP>0 , repeat sets 2E and F NOVP times; otherwise omit 2E> F. 
Note if NOVP> 1 , YVP values must be strictly increasing. 



5 41-50 


FI 0.0 


51-60 FI 0. 0 



FI 0.0 


11-20 110 



6F10.0 XVP(I) 
VP (I) 


Station at which velocity 
profile is specified. If 
NOVP=1, any value may be 
used (e.g. , 0) . 

Number of points defining 
velocity profile at YVP. 
(NVPP>2) 


NVPP pairs of values 

XVP = fraction of chord 
VP = velocity profile value 

XVP (I) <XVP (1+1) 

XVP (1) <0. 

XVP (NVPP) >1.0 


(3) FILE DEFINITION 


3 

RHOB 

NAMELIST 


Default for all filenames is 0 
File name in the form = n, 

(1 <n<5) 




MIFILE 

Modal input file name 




CMFILE 

C-matrix file name 
































COLUMNS 


VARIABLE ; 

DESCRIPTION 

DP FILE 

Delta pressure file name 

SGFFILE 

Sect. gen. force file name 

GFFILE 

Generalized force file name 
Default for all file initial 


positions is 1 

MIF1 

Initial file position of MIFILE 

CMF1 

Initial file position of CM FILE 

DPF1 

Initial file position of DPFILE 

SGFF1 

Initial file position of 


SGFFILE '{ 

GFFT 

Initial file position of GFFILE 

i 

Note that CMFILE must be dis- 
crete; the others may be dis- 


Crete or the same. If any of 
DPFILE, SGFFILE or GFFILE are 
the same, note that files are 
positioned in the order. 
DPFILE/DPF1 , SGFFILE/SGF1 , and 


GFFILE/GFF1 


(4) PRINTED OUTPUT DEFINITION 


4 

RHOC 

NAMELIST 


The default for all output 
control is 0. 

The input may be of form 
n<0 - print fcr all conditions 
=0 - no print-out 
n>0 - print for first n 
conditions 




DPPRT 

SGFPRT 

GFPRT 

Delta pressure print control 
Sectional gen. force print 
control 

Generalized force print 
control 














CARD 

SET 

FIELD NO. 

FORMAT 

VARIABLE 

DESCRIPTION 

COLUMNS 




DWMPRT 

CMPRT 

PCMPRT 

Downwash matrix print control 
C-matrix print control 
Pressure coefficient matrix : 
print control 


(5) PLANFORM DISTRIBUTIONS 


5a 

RHOD 

NAMELIST 

NDWC 

Number of downwash chords 
0< (NDWC | * | NPDWC | <72 
(No Default) 

If NDWC<0, a default cosine 
distribution will be 
generated. 

YDWC (I)=C0S (i 7T / (2NDWC+1) ) 
XPWD (I) ~-COS (2iir / (2NPDWC+ 1) ) 


■ i ■ ■■ ; 


NPDWC 

Number of points per downwash 
chord. 

0< | NPDWC | <8 (No default) 

If NPDWC<0 , the user input in 
set 5B is assumed to be in 
local non-dimensional 
coordinates (EAR notation) 



• : 

NPRC (1) 

Number of pressure output 
chords . ■ ■ :*'*"' * 

If NPRC <0, the user input in 
set 5C is assumed to be in 
local non-dimensional 
coordinates 




NPPRC 

Number of points per pressure 
output chord- 

(1) If no user output chords are specified, but printed output 
or file output is indicated, a default set of 11 output 
chords (and for pressure, 21 output points/chord) is used 

(U = .01,. 1, .2,.. .,.8,. 9,. 99) , (£ = -.99, -.9, -.8,.. .,.8, 

.9, .99) . 




NSGFC (1) 

. 

' ' J 

Number of sectional force 
output chords. If NSGFC<0, 
the user input in set 5D is 
assumed to be in non- 
dimensional coordinates. 
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CARD 

FIELD NO. 

FORMAT 

VARIABLE 

DESCRIPTION 

SET 

COLUMNS 





If NDWCCO, omit set 5B. The specified default cosine distributions 
will be used. Otherwise, repeat set 5B NDWC times. Note that input 
in 5B should be in physical coordinates if NPDWC>0, and in local non- 
dimensional coordinates if NPEWCCO . 


5B 

mm 

F10.0 

YDKC 

Y or Y of downwash chord 



6F10.0 

XDWP(I) 

X or X of points on downwash 





chord 


HU 

(7F10.0) 


- 1= 1 , NPDWC 


If NPRC=0, omit set 5C; otherwise, repeat set 5C NPRC times.. Note 
that input in 5C should be in physical coordinates if NPRC > 0, and 
in local non-dimensional coordinates if NPRC < 0. 


5C 

BHBQ 

F10.0 

YPRC ' “ . 

Y or _Y 

pressure output chord 



6F10.0 

XPPT (I) 

X or X 

points on pressure 

!' ' 


(7F10.0) 


output 

chord. 1=1, NPPRC 


If NSGFC=0, omit set 5D. Note that input in 5D should be in physical 
coordinates if NSGFC>0, and in local non-dimensional coordinates if 
NSGFC<0 . 


5D 

LIST 

7F10.0 

YSGFC (I) 

Y or Y of sectional force 





output chords. I=1,NSGFC 


(6) CONDITION DEFINITION 


6 

RHOE 

NAMELIST 

RVALUE 

List of reduced frequencies, 
{*<0) 




MACHN0 

List of Mach numbers 
(0<M<1.0) 

Both items must have at least 
one value, and may have up to 
20 values . 


(7) MODAL INPUT DEFINITION 


7A 

RHOF 

NAMELIST 

NIPTS (I) 

Number of modal input 

points/ 





input zone NIFTS(I)>3, 

1=1, 





NOCS+1. Alternatively 

for 
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CARD 

SET 

FIELD NO. 

FORMAT 

VARIABLE 

DESCRIPTION 

COLUMNS 

7A 

(cont 

. 

. 

) 

. | . 

5 

. 1 

NDMDS 

ZF 

CRI (I) 

DGCST 

GPGUST 

GPGREF 

NOCS#0, all input points may 
be input in a single block, 
in which case the program 
will determine which input 
points lie in which input zone. 

Number of displacement modes 

A user supplied multipli- 
cative factor which will be 
applied to modes input in 
card set 7C. (Default is 1.0) 

Control rotation input indi- 
cator 

n<0. Delta 3Z/9X will be speci- 
fied at required hinge points 
n=0 , no input 
n>0. Cubic coefficients of 
control rotation will be 
specified 

(Default is 0 for all I) 

1=1 ,NOCS 

Di screte gust option 
(Default = .F.) 

DGUST = . T. will cause a dis- 
crete gust mode to be appended 
to the set of displacement 
modes. 

Gradual penetration gust 
option (Def ault=. F. ) GPGUST=.T. 
will cause a gradual pene- 
tration gust mode to be 
appended to the set of dis- 
placement modes. 

Gradual penetration gust 
reference (zero phase) point. 
(Default = 0) 


If MI FILE *0, i.e., if modal input is to be read from MIFILE, or if 
NDM DS=0, omit sets 7B and C. Otherwise, repeat sets 7B and C once 
or once per each input zone, i.e., either 1 or NOCS+1 times. 



CARD 

SET 

FIELD NO. 
COLUMNS 

FORMAT. 

VARIABLE 

DESCRIPTION 

■ 

LIST 

6F10.0 

X(I) , Y (I) 

Modal input points, 1=1 , NIPTS 
note (X (I) ,Y(I) ) * (X (J) , Y ( J) ) 
for I#J . * . 

Repeat set 7C once for each displacement mode, i. e, NDMDS times. 

7C 

LIST 

7F10.0 

Z (I) 

Modal displacement at point 
<X(I) ,Y(I)) ,I = 1,NIPTS 

If any CRI (I) #=0 

, repeat set 7D once for each CRI (I) #0. 

7D 

RHOG 

NAMELIST 

• 

CCR (I, J) 
or 

DZDX(I,J) 

Cubic coefficients of control 
rotation (Default =0) 

1=1,4; 1-C 0 , 2-C lr 3-C z , 4-C 3 
J=1, NDMDS 

Only those terms which are 
non-zero need be input 

Delta DZ/DX at control hinge 
equation points, 1=1,4. 

(U = 0,1/3, 2/3,1) 

J= 1 , NDMDS 


(8) TERMINATION DEFINITION 


8 

1 

1-7 

A7 , 3 X 

LNAME 

Termination indicator 
[Blank 1 

["EXIT” J CALL EXIT 

"RETURN" execute return to 
calling overlay 






Anything else, execute 






CALL OVERLAY (LNAME, LI ,L2 , 0) 


2 

1 1-20 

110 

LI 

Primary level overlay no. 


3 

21-30 

1 10 

L2 

Secondary level overlay no. 

























3.8 PROGRAM OUTPUT 

3.8.1 Program Result s 

Printed output of program results consists of an initial block 
of information reflecting the user's input, followed by those 
intermediate and final results specified in card set 4, Printed 
Output Definition. The output controlled by card set 4 is 
calculated for each k-value Mach number condition; the user may 
elect to have all, none, or some first n conditions printed. 

The intermediate output consists of downwash matrices, C-matrices, 
and coefficients of assumed main surface pressure terms. If 
C-matrix printout is specified and control surfaces exist, the 
C-matrices for all controls as well as the main surface pressure 
terms are printed. If downwash matrix printout is specified 
and control surfaces exist, the residual downwash as well as 
full downwash will be printed. Note that the output of full 
and residual downwash is cyclic, four modes at a time. 

The final results consist of unsteady pressures, sectional and 
total generalized forces. Generalized force output consists 
of a single matrix per condition- Sectional force output consists 
of a matrix per sectional force chord (in the order specified) 
per condition. Unsteady pressure output consists of 
real/imaginary and amplitude/phase per output point (in the order 
specified) written two modes at a time for all pressure modes 
per condition. 

Following all condition output, a summary of the CMFILE index 
is given if CMFILE is present. Finally a summary of maximum 
core required and central processor time used is given. 

All normal output includes the user's run title with date appended 
and k-value and Mach number identified (if applicable). 

Binary output from the program consists of all unsteady pressure 
results if DPFILE^O, all sectional force results if SGFFILE^O, 
and all generalized force results if GFFILE/O. The fcrm of the 
information is described in section 3 .5 (Note that a user may 
have all binary output placed on the same file, in which case 
the results are interspersed on a condition basis! 

3.8.2 Program Diagn os tic s 

Program diagnostics may occur during input preparation or 
execution of the problem. The RHOIV input processor attempts 
to read and check all user input, identifying as many data errors 
as possible. If any errors are discovered during input processing 


the execution is terminated following input. Errors which may 
be discovered include: 

(1) exceedance of program size restrictions, 

(2) illegal planform definition, 

(3) illegal distribution of points on the planform.,*,..-. 

(4) inaccurate file specifications, 

(5) illegal k-values or Mach numbers, 

(6) insufficient modal input definition. 

With one specific exception, if a user's input is processed with 
no errors, and sufficient memory allocation is provided, the 
job should be completed without user origin errors. If any 
errors do occur they should be of program or system origin. 

The exception to the above can occur during C-matrix calculations 
If a downwash chord is placed too close to the planforir tip or 
a control surface side edge, the chordwise integration grid 
routine CGRID may generate an illegal sequence of integration 
regions. If the condition £ = 3. 76|Y o |/S<.001 occurs a warning 
message will be printed to indicate a possible problem may occur 
for a particular integration chord. The condition will occur 
anytime a downwash chord is less than -02S from the tip or a 
control side edge. 

The message does not indicate an error, only the fact that the 
above condition has not been met. 

The following is a list of input processor diagnostics; the input 
numbers used are examples, 

P LANF ORM DEFINITION ERROR: 

1. Illegal no. leading edge defn. pts. (2<NLE<1 0) ,NLE =11 

2. Illegal no. trailing edge defn. pts. (2<NTE<10) ,NTE = 11 

3. Illegal no. control surfaces (0<NOCS<6) ,NOCS =7 

4. Illegal no. velocity profiles (0<NOVP) ,NOVP = -1 

5. YLE (1 ) YTE (1) 

YLE(NLE) / YTE(NTE) 


6 . 


7. Illegal leading edge definition, YLE not strictly increasing 

8. Illegal trailing edge definition, YTE not strictly increasing 

9. Leading edge and trailing edge intersect 

10. Control surface 1 does not lie within defined planform 

11. Control surface 2 and 1 are incompatible 

12. Velocity profile stations are not strictly increasing 

13. Illegal no. of velocity profile points for station 1 
(2<NVPP) ,NVPP = -1 

14. Velocity profile points, XVP, for station 1 are not strictly 
increasing 

15. *** SCAMP4 *** error *** on velocity station 1 non user 


error 

F ILE DEFINITI ON ERROR : 



File spacing error 4, encountered for initial position of 
DPFILE 

CMFILE name is not descrete 


PLAN F ORM DISTR IB UTION ERROR: 

1. Illegal no. downwash chords, (1<ABS (NEWC*NPDWC) <72) , NDWC = -10 

2. Illegal no. points/downwash chord, (1 <ABS (NPEWC) <8) ,NPDWC = 10 

3. Downwash point 1, not on defined planform 

4. Downwash chord 1, coincides with a control surface side 

edge 

5. Downwash point 2, coincides with a control surface hinge 

6. Downwash points no. 2 and no. 1, are coincident 

7. Pressure report chord 1, not on defined planfcrm 

8. Pressure report point for chord no. 1, not on planform 

9. Pressure report point 1, coincides with a control hinge 
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10. Pressure report chord 1 , coincides with a control side edge 

11. Sectional force report chord 1 , rot on defined planfcrm 

12. Sectional force report chord 1, coincides with a control 
side edge 

CONDITION D EFINITION E RROR ; 

1. Illegal no. k- values, (0 < no.KVAL) 

2. Illegal no. Mach nos., (0 < no. MACH) 

3. Illegal k-value, (0 < KVAL) 

4. Illegal Mach no. , (0 < MACH < 1.0) 

MODA L INPUT DE FINITION ERROR: 

1. Illegal no. displacement modes, (0<NDMDS) ,NDMDS = -3 

2. Illegal no. input points, (3<NIPTS), NIPTS=2 

3. Insufficient input pts. in zone 1- (3<NPZONE) ,NPZCNE =2 

4. Input pts. no. 2 and no. 3, are coincident 

5. File spacing error 1, on MIFILE 

6. I/O error MROW<MROWS on MIFILE while reading X,Y for rone 1 

7. I/O error MROW<MROWD on MIFILE while reading Z for zone 1 

8. I/O error MROW<MROWD on MIFILE while reading IIA for zone 1 



3.9 SAMPLE INPUT/OUTPUT 



The following input data, figure 4, is for the lifting surface 
shown in figure 5, page 76. A portion of. the output generated 
by this data case is given on pages 77-100- 





RHOIV SAMPL E DATA INPUT 
CARD COLUMNS 


SET 

1 

11 

21 

31 

41 

51 

1 

NASA TM 

X-2909 leading/trailing edge 

controls 


2A 

$RHOA MSID=7000, 

NLE=2, NTE= 

2, NOCS= 

2/ NOVP-l $ 


2B 

0.0 

0.0 

1.539 

1.270 



2C 

1.763 

0.0 

1.763 

1.270 



2D 

LECS 

1.2165 

.923 

1.072 

1.0 



TECS 

1.628 

.923 

1.665 

1.072 


2E 

0.0 


41 




2F 

0.0 

.8 

.005 

.86818 

.011 

.95 


.015 

.965 

.02 

.97 

.025 

.9746 


.0375 

.9819 

.05 

.9 877 

.0625 

.992 


.9875 

1.0001 

.1125 

1.0065 

.1375 

1.0017 


.1625 

1.0161 

.1875 

1.0199 

. 2125 } 

1.0231 


.2625 

1.0281 

.3125 

1.0321 

.3625 f 

1.0349 


.4125 

1.0367 

.4625 

1.0377 

. . 5125 j\- 

1.0379 


.5625 

1.0373 

.6125 

1.0359 

.6625 - 

1.0336 


.7125 

1.0303 

.7625 

1.0259 

. 7875a 

1.0231 


.8125 

1.0199 

.8375 

1.0161 

.8625 

1.0117 


.8875 

1.0065 

.9125 

1.0001 

.9375 

.992 


.95 

.9877 

.9625 

.9819 

.975 

.9746 


.98 

.97 

.985 

.965 

.989 

.95 


.995 

.86818 

1 . 

.8 




3 $RHOB CMFILE=1 $ 

4 $RHOC DPPRT=-1,CMPRT=-1,DWMPRT=-1,PCMPRT=-1 ,SGFPRT=-1,GFPRT=-1 $ 

5A JRHOD NDViC=-9, NPDWC=7 $ 

6 $RHOE RVALUE (1) =0. , MACHNO (1 ) =0 . 8 $ 

7 A $RHOF NIPTS (1) 7, NDMDS=5 $ 

7B 0.0 0.0 1.212 1.0 1.2165 .923 

1.3971 1.071 2.628 .923 1.665 1.072 

1.763 1.0 

7C 1. 1. 1. 1. 1. 1.1. 

.8815 -.3305 -.335 -.5156 -.7465 -.7835 -.8815 

0. .098 0. 0. 0. 0. 0. 

0. 0. 0- 0. 0. 0. -.1158792 

0. .098 0. 0. 0. 0. -.1158792 

8 EXIT 


FIGURE 4 
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CALCULATE UNSTEADY LOADINGS CAUSED BY NOTIONS OF A LIFTING 
SUOFACE with leading and/or trailing eoge, sealeo gap 
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.8000 .8662 • .9500 .9650 .9700 .9766 .9819 .9877 .9920 1.0001 1.0065 1.0117 1.0161 1.0199 1.0231 1.0261 

1.0321 1.0369 1.0367 1.0377 1.0379 1.0373 1.0359 1.0336 1.0303 1.0259 1.0231 1.0199 1.0161 1.0117 1.0065 1.0001 

.9920 .9877 .9819 .9766 .9700 .9650 .9500 .6682 .8000 
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4.2 COMMON BLOCK -USAGE 

The RHOIV program uses both BLANK and LABELED common. 

The LABELED common blocks are used for communication between 
the primary and secondary overlays, and for communication between 
routines in a secondary overlay. The block names and contents 
are described on the following pages. 

The T heading on the following pages refers to variable type: 

I - Integer, R - Real, C - Complex, L - Logical, H - Hollerith. 

BLANK common is used in most secondary overlays as a variable 
length working area. In general the program of a overlay 
calculates the area required for arrays in the various subroutines 
and passes a dimension and first word address of each array 
through the subroutine calling sequence. A description of the 
area used by each overlay is given in section 3.3. 



LABELED COMMON NAME: BASIC LOAD LEVEL: PRIMARY (1,0) 

DESCRIPTION: Basic information, used by all modules. 


VARIAELE 

a 

DIM. 

ENG. NOM. 

DESCRIPTION 

INDCM 

i 


- 

C- MATRIX CALCULATION 
INDICATOR, 0-MAIN SURFACE, 
>0- CONTROL SURFACE 

SYM 

i 



ANALYSIS SYMMETRY INDICATOR 
1- SYMMETRIC 





2- ANTISYMMETRIC 

BO 

R 


b 0 

REDUCED FREQUENCY REFERENCE 
LENGTH 

S 

R 


s 

PLANFORM SEMI -SPAN 

YROOT 

R 



Y VALUE OF PLANFORM ROOT 
STATION 

kval 

R 


k = oj/V 

K-VA1UE, REDUCED FREQUENCY 

kvalr 

R 


k r = b 0 w/V 

REFERENCE K-VALUE 

KSQD 

R 


k 2 


MACH 

R 


M 

MACH NO. 

BETA 

R 


3 = 


BETAS QD 

R 


3 2 


RTITLE 

I 

9 


USER RUN TITLE, WITH DATE 
APPENDED 
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LABELED COMMON; NAME; CMLIB LOAD LEVEL; PRIMARY (1,0) 


DESCRIPTION: Variables describing the contents of CMFILE are 

maintained in CMLIB during the condition cycle. 


NO. 

VARIABLE 

T 

DIM. 

. ENG. NOM. 

DESCRIPTION 

1 

FINDCM 

| 


■ - 

INDICATOR TO SEARCH FOR A 
C-MATRIX 

2 

CMFOUND 

1 



INDICATOR FOR C-MATRIX FOUND 

3 

IDNAME 

I 



INDEX MATRIX NAME = 10HCMFIL 
INDEX ON WRITE, CHECKED ON 
READ 

4 

NOCM 

I 


•' | 

NUMBER OF C-MATRICES ON 
CMFILE 

5 

NMSNTRY 

I 



NUMBER OF MAIN SURFACE 
ENTRIES IN INDEX,. (=NO. MAIN 
SURFACE C-MATRICES ON CMFILE) 

6 

CDATE 

I 



LIERARY CREATION DATE 

7 

LMDATE 

I 



LAST DATE C-MATRICES ADDED 

8 

LINDEX 

I 



LENGTH OF CMFILE INDEX 

9 

MSENTRY 

I 



POSITION WITHIN INDEX OF A 
MAIN SURFACE ENTRY 

10 

NCSE 

I 



NUMBER OF CONTROL SURFACE 
ENTRIES ASSOCIATED WITH A 
MAIN SURFACE ENTRY 

11 

MATPOS 

I 



MATRIX POSITION WITHIN 
CMFILE . 

12 

DATE 

I 



CURRENT DATE 









LABELED COMMON 

NAME: CMVAL 

LOAD LEVEL: . SECONDARY ( 1 , 5) 

DESCRIPTION: 

Used by most C-Matrix routines. 

NO. 

VARIABLE 

T 

r — ~ — — — — — 
DIM. 

ENG. NOM. 

DESCRIPTION 

1 

ROWC 

C 

(72,8) 

. . , C 

ROWS OF C-MATRIX ASSOCIATED 
WITH DCWNWASH POINTS ON CHORD 

2 

GXYY 

C 

(72,8) 

G(x,y,y) 

SUBTRACTION TERMS ASSOCIATED 
WITH EVALUATION OF DIPOLE 
SINGULARITY 

3 

GPXYY 

C 

(72,8) 

G' (x,y,y) 


4 

CIPK 

c 

16 


CHORDWISE INTEGRAL OF PRESSURE 
TERMS TIMES KERNEL 

5 

KERN 

c 


K (x 0 ,y Q/ k,M) 

KERNEL FUNCTION VALUE 

























LABELED COMMON NAME: 


COUNT 


LOAD LEVEL: PRIMARY (1,0) 


DESCRIPTION: 


Miscellaneous counters used during condition cycle 
to size working storage. 


ENG. NOM. 


DESCRIPTION 



NSPT 


NCPT 


NPTRM 


NDMDS 


NPMDS 


NZONES 


MIPTS 


MIPSGF 


MICGF 


MIPGF 


NOVP 


LOCVP 



NUMBER OF SPANWISE MAIN 
SURFACE PRESSURE TERMS 

NUMBER OF CHORDWISE MAIN 
SURFACE PRESSURE TERMS 


TOTAL NUMBER OF ASSUMED 
MAIN . SURFACE PRESSURE TERMS 

NUMBER OF DISPLACEMENT MODES 

NUMBER OF PRESSURE MODES 
(NOTE, EXCEPT IN THE CASE OF 
A GUST ANALYSIS, N DM DS= NPMDS. 
FOR A GUST ANALYSIS 
NPMDS= NDMDS ♦ 1 . ) 

NUMBER OF MODAL INPUT ZONES 

MAXIMUM NUMBER OF INPUT 
POINTS IN ANY MODAL INPUT 
ZONE 

MAXIMUM NUMBER OF QUADRATURE ' 
POINTS REQUIRED FOR CHORDWISE 
INTEGRATION FOR ANY SECTIONAL 
FORCE OUTPUT CHORD. 

MAXIMUM NUMBER OF QUADRATURE 
CHORD REQUIRED FOR SPANWISE 
INTEGRATION FOR TOTAL GENERAL. 
FORCES. 

MAXIMUM NUMBER OF QUADRATURE 
POINTS REQUIRED FOR CHORDWISE 
INTEGRATION ALONG ANY 
QUADRATURE CHORD FOR TOTAL 
GENERALIZED FORCES. 

NUMBER OF VELOCITY PROFILES 

LOCATION OF VELOCITY PROFILE 
INFORMATION IN BLANK COMMON 
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LABELED COMMON 

NAME: 

CQUAD 

LOAD LEVEL: SECONDARY (1,5/7) 

DESCRIPTION: 

Used by CHDINT and Chordwise Quadrature routines 
in C-xnatrix calculation. DELXI is also used during 
sectional and total generalized force integration. 

NO. 

VARIABLE 

T 

DIM. 

ENG. NOM. 

DESCRIPTION 

1 

XILIL 

R 



; «L ' 

£ OF CHORDWISE INTEGRATION 
REGION LOWER LIMIT 


XIDIL 

R 




K OF CHORDWISE INTEGRATION 
REGION UPPER LIMIT 


DELXI 

R 



«U-«L 



XIMID 

R 



( V ^/2 




Ill 












' “ — ' 1 * ■ ■■ 1 • — — 11 ' ' ' "I 

LABELED COMMON NAME: CSVAL LOAD LEVEL: SECONDARY (1, 1-3/5/7) 


DESCRIPTION: 

Used 

most 

by TETERM, DTETERM, LETERM, DLETERM, 
terms generated by CSINIT. 

' • 

O 

53 

VARIABLE 

T 

DIM. 

ENG. NOM. 

DESCRIPTION 

1 

INDCS 

I 



.*'*.• 

INDICATOR FOR TYPE CONTROL 
SURFACE,.. 1-FULL, : 2-TIP, 

3 -MID, 4- PARTIAL 

2 

DECS 

L 

• 




INDICATOR FOR LEADING EDGE 
CONTROL, . T. =LEADING EDGE, 
. F.=TRAILING EDGE. 

3 

SYMSF 

R 



S f 

SYMMETRY SIGN FACTOR 
♦1.0 - SYMMETRIC ANALYSIS 
-1.0 - ANTISYMMETRIC ANALYSIS 

4 

YCSI 

R 




INBOARD CONTROL SIDE EDGE 

5 

XHCSI 

R 



X. 

1 

INEOARC CONTROL HINGE POINT 

6 

XLCSI 

R 



; / 

LEADING EDGE VALUE AT INBOARD 
CONTROL SIDE EDGE 

7 

YCSO 

R 




OUTBOARD CONTROL SIDE EDGE 

8 

XHCSO 

R 



x o 

OUTEOARD CONTROL HINGE POINT 

9 

XLCSO 

. 

R 




LEADING EDGE VALUE AT OUT- 
BOARD CONTROL SIDE EDGE 

10 

DELYCS 

R 



0 

1 

H- 

CONTROL SPAN 

11 

TANLH 

R 



tanA H 

TANGENT OF HINGE SWEEP ANGLE 

12 

• 

BETAH 

R 



S H 


13 

BHSQD 

R 



b h 


14 

TANLL 

R 



tanA L 

TANGENT OF LEADING EDGE SWEEP, 
FOR A LEADING EDGE CONTROL 

15 

BETAL 

R 





16 

ELSQD 

R 



6 L 


17 

XIE1 

R 




PHYSICAL, OR CONSTANT PERCENT 
CHORE EXTENSION OF HINGE FOR 









LABELED COMMON NAME: DWPTS LOAD LEVEL: PRIMARY (1,0) 


DESCRIPTION:: Used in INPREP, DWPREP, RDWRTC, and CMCALC 


NO. 

VARIABLE 

T 

DIM. 

ENG. NOM. 

DESCRIPTION 

1 

NDWC 

I 



NUMBER OF DOWNWASH CHORDS 

2 

YDWC 

R 

72 

. Y 

DOWNWASH CHORDS , S PANWI SE 
COORDINATE VALUES FOR DOWNWASH 
POINTS 

3 

NPDWC 

I 



NUMBER OF POINTS PER DOWNWASH - 
CHORD 

4 

NDWP 

I 



TOTAL NUMBER OF DOWNWASH 
POINTS (=NDWC*NPDWC) 

5 

XDWP 

R 

72 

x 

DOWNWASH POINT STREAMWISE 
COORDINATE VALUES 



LABELED COMMON NAME 


LOAD LEVEL : PRIMARY (1,0) 


DESCRIPTION: 


Read in INPREP - used by RHOTV. 
Specifies termination procedure. 


NO. VAR I A RLE 


ENG . NOM . 


DESCRIPTION 


LNAME 


2 Li 


3 L2 



TERMINATION INDICATOR 


( BLANK 1 
"EXIT" J 


CALL EXIT 


= "RETURN", EXECUTE RETURN 
OTHERWISE 

CALL OVERLAY (LNAME, LI, L2,0) 
PRIMARY LEVEL OVERLAY NO. 
SECONDARY LEVEL OVERLAY NO. 













LABELED COMMON 

NAME: FILES 

LOAD LEVEL: PRIMARY (1,0) 

DESCRIPTION: 

Defines all 

files used within the 1 program. . 

NO. 

VARIABLE 

T 

DIM. 

ENG. NOM. 

DESCRIPTION 

1 

MIFILE 

I 



MODAL INPUT FILE NAME 

2 

CMFILE 

I 



C- MATRIX FILE NAME 

3 

DPFILE 

I 



DELTA PRESSURE FILE NAME 

4 

SGFFILE 

I 


; ■ ■ . 

SECTIONAL FORCE FILE NAME 

5 

GFFILE 

I 



GENERALIZED FORCE FILE NAME 






NOTE 1-5 ARE USER ASSIGNED 






NAMES. 

6 

IN 

I 



STANDARD INPUT FILE NAME 

7 

OUT 

I 



STANDARD OUTPUT FILE NAME 

8 

RHOSC1 

I 



SCRATCH FILE NO. 1 NAME 


(MISFILE) 

I 



MODAL INFORMATION SCRATCH FILE 


(RES FILE) 

I 



NAME, RESULT SCRATCH FILE NAME 

9 

RHOSC2 

I 


•V - : ■■ 

SCRATCH FILE NO. 2 NAME 


(INSFILE) 

I 



INPUT SCRATCH FILE NA-IE 


(CMSFILE) 

I 



C-MATRIX SCRATCH FILE NAME 


(COFFILE) 

I 



COEFFICIENT FILE NAME 

10 

MIF1 

I 



INITIAL FILE POSITION OF MIFILE 

11 

CMF1 

I 



INITIAL FILE POSITION OF CMFILE 

12 

DPF1 

I 



INITIAL FILE POSITION OF DPFILE 

13 

SGFF1 

I 



INITIAL FILE POSITION OF SGFFIL 

14 

GFF1 

I 



INITIAL FILE POSITION OF GFFILE 



LABELED COMMON 

NAME: 

. FILES LOAD LEVEL: PRIMARY i 

(1,0) 

DESCRIPTION: 

Field 

length control/utilization: inf oriration. 




DIM. ENG. NOM. 


DESCRIPTION 


CURRENT PROGRAM FIELD LENGTH 

MAXIMUM PROGRAM FIELD LENGTH 
USED TO CURRENT TIME 

INITIAL PROGRAM FIELD LENGTH 

JOB CARD FIELD LENGTH (I.E., 
MAXIMUM ALLOWABLE FIELD LENGTH. 

MINIMUM FL REQUIRED FOR INPREP 

MINIMUM FL REQUIRED FOR MIPREP 

MINIMUM FL REQUIRED FOR RESPREE 

MINIMUM FL REQUIRED FOR RDWRTC 

MINIMUM FL REQUIRED FOR CMCALC 

MINIMUM FL REQUIRED FOR SOLN 

MINIMUM FL REQUIRED FOR RESULTS 

LENGTH OF INTERPOLATION 
INFORMATION ARRAYS 

LENGTH OF VELOCITY PROFILE 
INFORMATION 

LENGTH OF CONTROL ROTATION 
COEFFICIENTS INFORMATION 














LABELED COMMON NAME: KRNTERM 


DESCRIPTION: Used by kernel function 


NO. 

VARIABLE 

T 

: DIM. . 

' ENG . NOM. . 

1 

FK 

R 


fk 

2 

FKSQD 

R 


(fk) 2 ‘ 

3 

RFKSKYS 

R 


/(fk) i +(ky o ) 

4 

HK 

R 


hk 

5 

HKSQD 

R 


(hk) 2 

6 

RHKSKYS 

R 


/(hk) 2 + (ky Q ) 

7 

HOYO 

R 


h /y 0 

8 

aproxr 

R 



9 

APFOXI 

R 



10 

SIR 

R 



1 1 

S1I 

R 



12 

S2R 

R 



13 

S2I 

R 



14 

S3R 

R 



15 

S3I 

R 
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LABELED COMMON NAME 


KRNVAR 


LOAD LEVEL: SECONDARY (1 , 5) 


DESCRIPTION: Used by C-matrix chcrdwise quadrature routines, and 

kernel function. 


VARIABLE 


ENG. NOM.. 


IDWC 

' I 

Y 

R 

IPDWC 

I 

IDWP 

I 

x 

R 

YO 

R 

YOSQD 

R 

KYO 

R 

KYOSQD 

R 

B YOSQD 

R 

EK YOSQD 

R 

XO 

R 

KXO 

R 

CMACH 

R 


DESCRIPTION ' 


DOWNWASH CHORD NO. 

Y VALUE OF DOWNWASH POINT 
POINT/DOWNWASH CHORD 
DOWNWASH POINT NO. 
DOWNWASH POINT 


(ky 0 ) 2 
(6y 0 ) 2 

(8ky n ) 2 


In (2-2M) /2 


CONSTANT FOR KERNEL CALCU- 
LATION 














LABELED COMMON NAME: LCSTERM 

LOAD LEVEL: SECONDARY (1,5/7) 

DESCRIPTION: 

Miscellaneous terms associated with full chord 
control pressure expressions. 

NO. 

' 

VARIABLE 

T 

DIM. 

ENG. NOM. 

DESCRIPTION 

1 

CIS 

C 

4 

C IS 

COEFFICIENT/SIDE EDGE FOR 
INVERSE SQUARE ROOT TERMS 

2 

CSI 

C 

4 

C S1 

COEFFICIENTS/SI EE EDGE FOR 
SQUARE ROOT TERM ASSOCIATION 
WITH FIRST SOLUTION 

3 

CS2 

c 

4 

C S2 

■' . 

SAME AS Cg! EXCEPT FOR SECOND 
SOLUTION 

4 

CL 1 

c 

4 

C L1 

COEFFICIENTS/SIDE EDGE FOR 
FIRST SOLUTION 

5 

CL2 

c 

4 

C L2 

SAME AS C L2 EXCEPT FOR SECOND 
SOLUTION 

6 

CAT 

c 

4 

C AT 

COEFFICIENTS/SIDE EDGE FOR 
ARC TANGENT TERM 

7 

GIS 

R 


G IS 

INVERSE SQUARE ROOT TERM 

8 

GS1 

. 

R 


G S1 

FIRST SOLUTION SQUARE ROOT 
TERM 

9 

GS2 

c 


g S2 

SECOND SOLUTION SQUARE ROOT 
TERM 

10 

GL1 

R 


g L1 

FIRST SOLUTION LOG TERM 

11 

GL2 

C 


g L2 

SECOND SOLUTION LOG TERM 

12 

GAT 

R 


g at 

ARC TANGENT TERM 

13 

CONI 

R 



MISCELLANEOUS CONSTANT TERMS 

14 

CON 2 

R 




15 

CON3 

R 




16 

CCON1 

C 




17 

CCON2 

C 












LABELED COMMON 

NAME: OPTIONS 

LOAD LEVEL: PRIMARY (1,0) 

DESCRIPTION: 

Miscellaneous options. 


NO. 

VARIABLE 

T 

DIM. 

ENG. NOM. 

DESCRIPTION 

1 

VELPFL 

L 

' 



VELOCITY PROFILE OPTION 
•T. = PROFILES EXIST 

2 

DGUST 

L 



DESCPETE GUST OPTION 

3 

GPGUST 

L 



GRADUAL PENETRATION GUST 
OPTION 

4 

GPREF 

R 



GRADUAL PENETRATION GUST 
REFERENCE POINT 

5 

IIAIN 

L 



INDICATOR FOR DIRECT INPUT OF 
II A PER MODAL INPUT ZONE 

6 

RESULT 

L 



INDICATOR FOR SOME RESULTS 
REQUIRED, I. E. EITHER UNSTEADY 
PRESSURE, SECTIONAL FORCES 
OR GENERALIZED FORCES. 

7 

SOLUTION 

L 



INDICATOR FOR SOLUTION 
REQUIRED, I.E. EITHER RESULT 
OR C- MATRIX, DOWN WASH MATRIX, 
ON PRESSURE COEFFICIENT MATRIX 
OUTPUT REQUESTED. 




LABELED COMMON NAME: 

PLNGEO 

LOAD LEVEL: PRIMARY (1,0) 

DESCRIPTION: 

Defines 

planform geometry. 

NO. 

r— — — 

VARIABLE 

T 

DIM 

. 

ENG. NOM. 

DESCRIPTION . .. 

1 

MS ID 

I 




MAIN SURFACE C-MATRIX I. D. 

2 

NLE 

I 




NUMBER OF LEADING EDGE DEFN. 
POINTS 

3 

XLE 

R 

10 



LEADING EDGE DEFN. PTS. 

. (STRAIGHT LINE SEGMENTS) 

4 

YLE 

R 

10 



5 

DXLEDY 

R 

9 



SLOPE OF LEADING EDGE SEGMENTS 

6 

NTE 

I 




NUMBER OF TRAILING EDGE DEFN. I 

7 

XTE 

R 

10 




8 

YTE 

R 

10 


' • •* ... • '• ' ' ' 

TRAILING EDGE DEFN. PTS. 
(STRAIGHT LINE SEGMENTS) 

9 

DXTEDY 

R 

9 



SLOPE CF TRAILING EDGE SEGMENTS 

10 

NOCS 

I 




NUMBER OF CONTROL SURFACES 

11 

XHLI 

R 

6 



CONTROL SURFACE HINGE 
INBOARD DEFIN. POINT 

12 

YHLI 

R 

6 




13 

XHLO 

R 

6 



CONTROL SURFACE HINGE 
OUTBOARD DEFN. POINT 

14 

YHLO 

R 

6 




15 

CXHLDY 

R 

5 



SLOPE OF CONTROL HINGE 

16 

CSID 

I 

6 



CONTROL SURFACE C-MATRIX I. D. 

17 

CSTYPE 

I 

6 



CONTROL SURFACE TYPE 

18 

CSRS 

I 

6 



CONTROL SURFACE RELATED 
SURFACE 

19 

CSAO 

I 

6 



CONTROL SURFACE AREA ORDER 

20 

CSRI 

I 

6 



CONTROL SURFACE INPUT INDICATOJ 

21 

XHLBI 

R 

6 



HINGE INBOARD SIDE EDGE DEFN. 
POINT VALUE IN BAR NOTATION 
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22 

XHLBO 

R 

6 


HINGE OUTBOARD SIDE EDGE 
DEFN. POINT VALUE IN BAR 
NOTATION 

23 

XLEI 

R 

6 


LEADING EDGE VALUE AT THE 
CONTROL INEOARD SICE EDGE 

24 

XLEO 

R 

6 


LEADING EDGE VALUE AT THE 
CONTROL OUTBOARD SIDE EDGE 




LABELED COMMON NAME: PRSPVAL 

LOAD LEVEL: SECONDARY (1, 1-3/5/ 

DESCRIPTION: 

Used to specify (£,n) values to ETACVAL and other 



routines. 

Describes current chord, control surface 



and pressure type. 


NO- 

VARIABLE 

T 

DIM. 

ENG. NOM. 

DESCRIPTION 

1 

ETA 

R 


n 

SPANWISE LOCATION OR COORDINATI 

2 

AETA 

"R 


In| 


3 

ETAB 

R 


n 

NON-DIMENSIONAL SPANWISE 






COORDINATE 

4 

XI 

R 



CHORDWISE (STREAMWISE) LOCATION 
OR COORDINATE. 

5 

XIB 

R 


. . c : 

NON-DIMENSIONAL CHORDWISE 






COORDINATE 

6 

DXIBE 

R 


8i/3n 


7 

ICS 

I 



CONTROL SURFACE NUMBER 
(=0 WHEN NO CONTROL SURFACE IS 
CONCERNED) 

8 

CPT 

L 



CONTROL PRESSURE INDICATOR, : . 
CPT=.T. IF CONTROL SURFACE 
PRESSURE TERM IS BEING 
CALCULATED. 

9 

XIL 

R 


£l<n> 

PLANFORM LEADING EDGE VALUE 
AT n 

10 

DXILE 

R 




11 

XIM 

. 

R 



PLANFORM MIDCHORD VALUE AT n 

12 

XIC 

R 


(n) 

PLANFORM CONTROL NO- ICS HINGE 





‘ - c - 

VALUE AT f| . (NOTE - THIS MAY 
BE A LINEAR EXTENSION.) 

13 

DXICE 

R 


3C c (n)/3n 


14 

XIT 

R 


C t (n) 

PLANFORM TRAILING EDGE 






VALUE AT n 

15 

DXITE 

R 


3C t (n)/9n 


16 

B 

R 


b(n) 

PLANFORM SEMI -CHORD VAlUE 
AT U 

17 

KIND 

I 



KERNEL INDICATOR, USED ONLY 
DURING C-MATRIX CALC. 

18 

CQTYPE 

I 



CHORDWISE QUADRATURE TYPE, 
USED DURING C-MATPIX CALCU- 
LATION, AND SECTIONAL AND 
TOTAL GENERALIZED FORCES. 
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LABELED COMMON NAME: 


LOAD LEVEL: S ECONDARY ( 1 , 3/5) 


DESCRIPTION: Used by PRSPREP, SGFPREP, GFPREP, and irost 

C-matrix calculation routines. 


NO. VARIABLE 


ENG. 


DESCRIPTION 


FETA 

R 

. 72 

f (T1) 

SPANWISE MAIN SURFACE 
PRESSURE TERMS 

FPETA 

R 

72 ' ; 

f,( n) "v 


NCPTERM 

I 


g(C»n) or 
3g (S f n) /3 rt 

NUMBER OF CHORDWISE PRESSURE 
TERMS 

CPTERM 

R 


CHORDWISE PRESSURE. TERMS 
















LABELED COMMON NAME: PRTCTL 

LOAD LEVEL: PRIMARY (1,0) 

DESCRIPTION: 

Controls all printed output in result and solution 
sections- 

NO- 

VARIAELE 

B 

DIM. ' 

ENG. NOM. 

DESCRIPTION ■ 

1 

DPPRT 

i 



DELTA PRESSURE PRINT CONTROL 


SGFPRT 

i 



SECT. GEN. FORCE PRINT CONT.' 


GFPRT 

i 



GEN. FORCE PRINT CONTROL 

D 

CMPRT 

i 



C-MATRIX PRINT CONTROL 

5 

DWMPRT 

i 



DOWNWASH MATRIX PRINT CONTROL 

6 

PCMPRT 

i 



PRESSURE COEFF. MATRIX PRINT 
CONTROL 

n<0 all cond. (on input) 

= 0 - no print 

n>0 - print first n cond. 
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LABELED COMMON 

NAME: QUADWTS 

LOAD LEVEL,: PRIMARY (1 ,0) 

DESCRIPTION: 

Used by GFSGP.ID, GFCGRID 

in prep, routines, and by 



CMCALC and chordwise quadrature routines in C-iratrix 



calculation 

• • 


NO. 

VARIABLE 

T 

DIM. 

ENG. NCM. 

DESCRIPTION 

1 

ALEG4 

R 

2 



2 

HLEG4 

R 

2 


PREFIX , 

3 

ALEG8 

R 

4 


A - ABSCISSAE 

4 

HLEG8 

R 

4 


H - WEIGHTS 

5 

AISQR5 

R 

5 


ROOT - TYPE QUADRATURE WEIGHT 
FUNCTION 

6 

HISQR5 

R 

5 


LEG - LEGENDRE 

7 

AISQR10 

R 

10 


ISQR - INVERSE SQUARE ROOT 
SQR - SQUARE ROOT 

8 

HISQR10 

R 

TO 


LOG - LOG 

9 

ALOG4 

R 

4 



10 

HLOG4 

R 

4 


4THR - FOURTH ROOT 

11 

ALOG8 

R 

8 


SUFFIX NO. CUADRATURE POINTS 

12 

HLOG8 

R 

8 



13 

ASQR5 

R 

5 



14 

HSQR5 

R 

5 



15 

ASQR10 

R 

10 



16 

HSQR 1 0 

R 

10 


# 

17 

A4THR5 

R 

5 



18 

H4THR5 

R 

5 



19 

A4THR10 

R 

10 



20 

H4THR10 

R 

10 
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LABELED COMMON 

NAME: SQUAD 

LOAD LEVEL: SECONDARY (1 ,5) 

DESCRIPTION: 

Used by CMCALC, SPNINT 


NO. 

VARIABLE 

B 

DIM. 

ENG. NOM. 

DESCRIPTION 

1 

SQTYPE 

i 



SPANWISE QUADRATURE TYPE 

2 

SQWT 

R 



SPANWISE QUADRATURE WEIGHT 

3 

ETALIL 

R 


n L 

n OF LOWER LIMIT FOR A SPAN- 
WISE INTEGRATION REGION 

4 

ETAUIL 

■ ' 

R 


n u 

n OF UPPER LIMIT FOR A SPAN- 
WISE INTEGRATION REGION 

5 

DELETA 

R 


n U' r, L 

SPANWISE INTEGRATION REGION 
WIDTH - 

6 

ETAMID 

R 


(T1 U +T1 L )/2 

REGION MIDPOINT 

7 

ALNDELY 

R 


Ln[.-(n 0 .-n L )/s-J 


8 

ETAL 

. ■ 

R 


n l 

D VALUE OF LOG SINGULARITY 
LOCATION DURING SPANWISE 
INTEGRATION. IT IS EITHER 
DOWNWASH CHORD, OR A 
CONTROL SIDE EDGE. 
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LAB EL ED G@J©lON 

NAME: TIMES 


LOAD LEVEL: PRIMARY (1,0) 

DESCRIPTION: 

Accumulated Cp times and 
various secondary level 

number of entries for 
overlays. 

NOi 

VARIABLE 

T 

DIM. 

ENG. 

NOM. 

DESCRIPTION 

1 

TINP 

R 




ACCUMULATED CP TIME IN INPREP 

2 

TMIP 

R 




ACCUMULATED CP TIME IN MIPREP 

3 

TRESP 

R 




ACCUMULATED CP TIME IN RESPREP 

4 

NTRWC 

I 




NO . ENTRIES INTO RDWRTC 

5 

TRWC 

R 




ACCUMULATED CP TIME IN RDWRTC 

6 

. ' 

NTCMC 

I 

5 



1) ASSOCIATED WITH MAIN 
SURFACE, k = 0 

2) ASSOCIATED WITH MAIN 
SURFACE, k > 0 

3) ASSOCIATED WITH CONTROL 
SURFACE, k = 0 

4) ASSOCIATED WITH CONTROL 
SURFACE, k > 0 

5) TOTAL 1) - 4) 

7 

. 

' 

TCMC 

R 

5 



ACCUMULATED CP TIME IN CMCALC 
(SAME AS NTCMC) 

8 

NTSOLN 

I 

- | 



NO. ENTRIES INTO SOLN 

9 

TSOLN 

R 

t 



ACCUMULATED CP TIME IN SOLN 

10 

NTRES 

I 




NO. ENTRIES INTO RESULTS 

11 

TRES 

R 




ACCUMULATED CP TIME IN RESULTS 
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4.3 PROGRAM AND SUBROUTINE DESCRIPTION 

Short abstracts for the various programs and subroutines in RHOIV 
are included in this section. A full description of each routine 
may be found in the program listing.' 

Routines are ordered alphabetically; their positional relation 
is shewn in figure 7, page 132 • 


PROGRAM - SUBROUTINE STRUCTURE 


BLOCK 

LENGTH 

BLOCK 

LENGTH 

BLOCK 

LENGTH 

(00,00) 

1 

(01,01) 

1 

(01,04) 

1 

ZRHOIV 

12720 

INPREP 

2 33 6 

RDWRTC 

1526 



ERRPRT 

1401 





PLNINCK 

320 



(01,00) 

1 

CSINCK 

607 

(01,05) 

- 1 

/BASIC/ 

24 

VP INCH 

167 

/CMVAL/ 

664 2 

/CMLIB/ 

14 

VPINCK2 

214 

/KRNVAR/ 

16 

/COND/ 

52 

SCAMP 4 

250 

/SQUAD/ 

10 

/COUNT/ 

20 

COMCU 

2512 

CMCALC 

620 

/DWPTS/ 

124 

CU3IC2 

120 

SGRID 

414 

/ENDIT/ 

3 

DERIV1 

107 

ORDER 

50 

/FILES/ 

16 

DERIV2 

103 

GXYYCAL 

432 

/PLNGEO/ 

222 

NAMBLD 

:• 45 

SPININT 

300 

/PRTCTL/ 

6 

DWINCK 

573 

CHDINT ■ 

552 

/OPTIONS/ 

7 

PRSINCK 

503 

CGRID 

1457 

RHOIV 

212 

SGFINCK 

256 

CQLEG4 

75 

IIARDR 

47 

MIINCK 

1704 

CQLEG8 

75 

/QUADWTS/ 

176 



CQISR5 

64 

BLKDATA 

3 

(01,02) 

1 

CQISR10 

64 

/CSVAL/ 

26 

MIPREP 

273 

CQLOG4 

141 

/PRSPVAL/ 

22 

IIACAL 

213 

CQLOG8 v; 

141 

ETACVAL 

146 

PLAT El 

434 

CQSQR5 

66 

AZONE 

151 

PLATET 

155 

CQSR10 

66 

AINTL 

166 

PLATEA 

205 

CQ4R5 

66 

PLATEO 

544 

CRCOEFF 

672 

CQ4R1 0 

66 

ZEROCOL 

65 

CPRDR 

1 03 

CPFCTXI 

6 1 

/PRSTERM/ 

43 

VPCALC 

144 

DTETERM 

1271 

SPFCTE 

174 

CTHETA 

74 

DLETERM 

1110 

CMSTERM 

114 



DLFCT 

135 

/CQUAD/ 

4 

i ‘ 


/KRNTERM/ 

17 

/CSTERM/ 

102 

* (01,03) 

1 

KRNFCT 

375 

TETERM 

1206 

RESPREP 

722 

BE SK I 

142 

E2 

30 

CCRRDR 

54 

BI1ML1 

167 

DE2 

101 

DWPREP 

714 

APROX 

232 

/LCSTERM/ 

140 

VPRDR 

103 

SERIES1 

115 

CSINIT 

422 

VPCALC 

144 

SERIES2 

170 

LETERM 

400 

PRSPREP 

440 

SERIES3 

142 

LFCT 

52 

SGFPREP 

677 

SICI 

216 

MATIO 

340 

GFPREP ' 

731 





GFSGRID 

307 





GFCGRID 

720 





(01,06) 

1 

(01,07) 

1 



SOLN 

214 

RESULTS 

411 



PCMDWN 

1417 

FORMDP 

1250 



SOLVE 

531 

RTHETA 

102 



PCMSPT 

323 

FORMQS 

772 





FORMQ ; 

725 


Figure 7 




( 1) SUBROUTINE AINTL (X, Y, NPTS , Z f NRCWZ , NCOL 1 ,NCOLS, SA, INDD, 

DZ1, DZ2) 

Interpolation cover routine for generation of displacements and 
slopes at unsteady aerodynamics control points - specified in 
the local axis system. 

(2) SUBROUTINE APROX 

Routine APROX performs the integration of 

TAU*EXP (I*KYO*TAU) VSQRT (l. + TAU**2) . In TAU from -~tO H. The 
function TAU/SQRT ( T. + TAU 5 * * 2 ) is approximated by a series of 
exponential, NASA technical report R-48, page 8, 

TAU/SQRT ( 1 . +TAU**2) = 1. +C1 *EXP (E1*TAU) + C2*EXP (E2*TAU) + 
C3*EXP (E3*TAU) *SIN (PI*TAU) 

(3) SUBROUTINE AZONE (Y, X, NPTS, INB, ZONE) 

Assign a modal input zone number to points on chord Y. 


(4) FUNCTION EESKl (X) 

BESK1 = Kl (X) , modified Bessel function of the second kind of 
order 1. Kl (X) is calculated with two polynomial expansions, 
one for X.LT.2, and one for X.GE.2. Equations taken from Handbook 
of Mathematical Functions, National Bureau of Standards, 1967. 
Bessel Functions of Integer Order , 9. 8 Polynomial Approximation, 
p. 378, 9.6.11 Ascending Series, p. 375. 


(5) FUNCTION BI1ML1 (X) 

BI1ML1 = (II (X) - LI (X) ) where II = Modified Bessel function 
of the first kind of order 1, LI = Struve function. 

Series expansion for I1-L1 is used for X.LE. 12.8, an asymptotic 
expansion is used for X.GT. 12.8. 

Equations taken from Handbook of Mathematical Functions, National 
Bureau of Standards, 1967. 

Struve Functions and Related Functions, 12.2.6 Asymptotic 
Expansion p. 498, 12.2.1 Power Series p. 498. 

Bessel Functions of Integer Order, 9.6.10 Ascending Series, P. 
378. 
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BLKDATA 



Defines all quadrature abscissae and weights. 



(7) SUBROUTINE CCRRDR (CCR, NPMOS) 

Read cubic coefficients of control rotation. 


(8) SUBROUTINE CGRID (NIR, XIQ, IRTYPE) 

Routine CGRID defines the chordwise integration schemes used 
to integrate chordwise pressure terms times one of the kernel 
expressions on a specified Chord. 

(9) SUBROUTINE CHDINT 

Perform the chordwise integration of chordwise pressure terms 
times kernel function. 


(10) OVERLAY (RHOIV, 1,5) PROGRAM CMCALC 


Calculate a C-matrix associated with main surface or control 
surface pressure terms for a set of downwash points at a 
particular k-value Mach number condition. 


(11) SUBROUTINE CMSTERM 

Calculate NCPTEFM ' chordwise pressure terms associated with the 
assumed main surface pressure terms. 


(12) SUBROUTINE CPFCTXI 

Calculate NCPTERM chordwise pressure terms 

KIND = GXYY, calculate G (XI, Y) , DG (XI , Y) /DETA 

= non-sing or singular calculate G (XI, ETA) 

CPT = .F.-raain surface terms, .T. -control surface terms 

LECS = . FALSE. -calculate only T.E. control surface terms 

.TRUE. -calculate both T.E. and L.E. control surface terms 
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(13) SUBROUTINE CQISR5 

Perform A 5 pt. inverse square root quadrature on chordwise 
pressure term times kernel function for XI over XILIL to XIUIL 
where the square root singularity is at XILIL. 

(14) SUBROUTINE CQISR10 

Perform a 10 pt. inverse square root quadrature on chordwise 
pressure term times kernel function for XI over XILIL to XIUIL 
where the square root singularity is at XILIL. 


(15) SUBROUTINE CQLEG4 

Perform a 4 pt. Gauss -Legendre quadrature on chordwise pressure 
term times kernel function for XI over XILIL to XIUIL. 


(16) SUBROUTINE CQLEG8 

Perform a 8 pt. Gauss -Legendre quadrature on chordwise pressure 
term times kernel function for XI over XILIL to XIUIL. 


(17) SUBROUTINE CQLOG4 

Perform a combination 4 pt. log and 4 pt. Legendre quadrature 
on chordwise pressure term times kernel function for XI over 
XILIL to XIUIL where the log singularity is at XILIL. 

i ' : . • • . . . . 

I- ■ • . ; ; . " ■ 

(18) SUBROUTINE CQLOG8 

Perform a combination 8 pt. log and 8 pt. Legendre quadrature 
on chordwise pressure term times kernel function for XI over 
XILIL to XIUIL where the log singularity is at XILIL. 


(19) SUBROUTINE CQSQR5 

Perform a 5 pt. square root quadrature on chordwise pressure 
term times kernel function for XI over XILIL to XIUIL where the 
square root slope singularity is at XILIL. 

NOTE: Abscissae and weights developed from Gauss-Nehler form, 

refer to Kopal, Numerical Analysis, pp. 381-2, using 
BCS program MEHQAH (Redman 1973) with ALPHA=0 , BETA= . 5 . 



(20) SUBROUTINE CQSQR10 

Perform a 10 pt. square root quadrature on chordwise pressure 
term times kernel function for XI over XILIL to XIUIL where the 
square root slope singularity is at XILIL. 

NOTE: Abscissae and weights developed from Gauss-Mehler form, 

refer to Kopal, Numerical Analysis, pp. 381-2, using 
BCS program MEHQAH (Redman 1973) with ALPHA=0 , BETA=. 5 . 


(21) SUBROUTINE CQ4R5 

Perform a 5 pt. fourth root quadrature on chordwise pressure 
term times kernel function for XI over XILIL to XIUIL where the 
3/4 root slope singularity is at XILIL. 

NOTE: Abscissae and weights developed from Gauss-Mehler form, 

refer to Kopal, Numerical Analysis, pp. 381-2, using 
BCS program MEHQAH (Redman 1973) with ALPHA=0 , BETA=. 25 


(22) SUBROUTINE CQ4R 10 

Perform a 10 pt. fourth root quadrature on chordwise pressure 
term times kernel function for XI over XILIL to XIUIL where the 
3/4 root slope singularity is at XILIL. 

NOTE: Abscissae and weights developed from Gauss-Mehler form, 

refer to Kopal, Numerical Analysis, pp. 381-2, using 
BCS program MEHQAH (Redman 1973) with ALPHA=0 , BETA=.25 


(23) SUBROUTINE CRCO EE F (NZONES, NDMDS, NPMDS, DELDZDX, CCR, 

1 Z, IIA, IIAP) 

Calculate (or read) the cubic coefficients of control surface 
rotation for subsequent use in the preparation routines. The 
coefficients are defined such that 

THETA (ETA) = CO + C1*EBCS + C2*EBCS**2 + C3*EBCS**3 
= change in streamwise modal slope, 

DELTA DZ/DX, at the hinge. 

EBCS = fraction (0.,1.) of control span 
(0. inboard, 1.0 outboard) 


(24) SUBROUTINE C THETA (DELDZDX, NDMDS) 



Given four equally spaced points on a control hinge, and DELTA 
DZ/DX at these points, determine the coefficients of the cubic 
which- will define DELTA DZ/DX, THETA, on the hinge. . 

(25) SUBROUTINE CSINCK (CSAREA, PLNERR, CSEPR) 

Read control surface input. Check legality. Assign control, type 
and determine related surfaces. 


(26) SUBROUTINE CSINIT 

Initialize variables for a control surface and generate 
coefficients of control pressure expression. 

(27) SUBROUTINE CLETERM (DGXIETA) 

Calculate the derivative of the pressure term G (XI, ETA) with 
respect to ETA for a leading edge control surface. 


(28) SUBROUTINE CTETERM (DGR, DGI) 

Calculate the derivative of the pressure term for a trailing 
edge control, DG (XI, ETA) /DETA = (DGR, DGI) 


(29) SUBROUTINE DWINCK (PLNERR, DWERR) 

Read (or use default definition for) downwash chords and points 
convert from bar notation if required, and check legality. 




SUBROUTINE EWPREP (NDWC, NPDWC, NDWP, NDWP2, NOCS, N ZONES , N DM DS , 
1 NPMDS, YDWC, XDWP , X,Y, ZONE, IFOS, NPZONE, 

Z,ZT, DZDX ,.DZDXT , W, WRI, IIA, IIAP, OCR) 


Form the basic downwash matrix, [W] = [DZ/DX + I*Z ], and write 
to RESFILE. The real portion of W is modified by any user 
supplied velocity profile, and a gust mode of the desired form 
is appended if requested. Any existing cubic coefficients of 
control rotation are copied from memory to RESFILE. 


(31) SUBROUTINE ERRPRT (IERRNO, ERRFLG, II, 12) 

Set error flag and print appropriate error message. 


(32) SUBROUTINE ETACVAL 

Given a spanwise station, ETA, determine the leading edge, 
trailing edge, and hinge line (if applicable) intersects, and 
calculate semi-chord and mid-chord values. 


(33) SUBROUTINE FORMDP (NPRC, NPPRC, NFPT, 

1 YPRC , XPPT, XBAR, 

2 CMSPT, FETA, AMPPHAS) 


NPTRM, 

PRESS, 


NOCS, NPMDS, 
MSPTRM, GXIETA, 


Calculate unsteady delta pressure at an indicated set of points 
(output points) for a k- value/Mach no. condition. 


DELP (X,Y ; J) = SUMf MSPTRM (X,Y; I) *CMSPT (I, J) , 1=1, NPTRM ] + 
j SUM[ FETA (Y,J) *GXIETA (X,Y;N) , N=1 ,NOCS ] 


(34) SUBROUTINE FORMQ (NDMDS, NPMDS, NPTRM, NOCS, MIC HD , MPICHD, 

1 Q, QMSPT, CMSPT, YICHD, GXIETA, XIPT, QTYPE, 

2 FETA, WZ) 

Calculate the generalized unsteady aerodynamic coefficient matrix 
(generalized forces) for a k-value/Mach number condition. 

Q (I, J) = INTEGRAL [ Z (XI ,ETA; I) *DELP (XI , ETA; J) *D[ XI ]*D[ ETA ] ] 

XI = XIL(ETA), XIT(ETA) 

ETA = 0, S 


(35) SUBROUTINE FORMQ S (NSGFC, NDMDS, NPMDS, NPTRM, NOCS, MPSGFC, 

1 Y, QS, QSMSPT , CMSPT, GXIETA,. XIPT, QTYPE, 

2 FETA, WZ) 
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Calculate the sectional generalized unsteady aerodynairic 
coefficient matrices (sectional forces) at an indicated set of 
chords for a k-value/Mach no. condition. 

QS(Y;I,J) = INTEGRAL[ Z (XI , Y; I) * (DELF (XI, Y; J) *D[ XI ] ] 

XI = XIL(Y) , XIT(Y) 

(36) SUBROUTINE GFCGRID (NIPTS, NPZONE, XI PT, CQWT, QTYPE) 

Determine the quadrature points: and associated weights to be 
used in integrating delta pressure times displacements along 
a chord for sectional or total generalized forces. 


(37) SUBROUTINE GFSGRID (NIC HD, YICHD , SQWT) 

Determine the quadrature chords and associated weights to be 
used in the s panwise integration cf delta pressure times 
displacements for total generalized forces. 

(38) SUBROUTINE GFPREP (MICHD, MPCHD, NZONES , NDMBS, NPMDS,. NSPT, 

1 NCPT, NPTRM, YICHD, SQWT, NIPTSZ, XIPT>- 

2 YIPT, CQWT, QTYPE, QMSPT, Z, WZ, G FETA, IIA, 

3 HAP, CCR, MICGF, MIPGF) 

Prepare information, independent of k-value and Mach no. , which 
will be used for calculation of generalized unsteady aerodynamic 
coefficients (generalized forces). The information is written 
to RESFILE. 


(39) SUBROUTINE GXYYCAL 

Calculate G(X,Y,Y) and GP(X,Y,Y) for all points on a downwash 
chord and initialize the associated C -matrix rows. 

(4) SUBROUTINE II ACAL (NZONES, MIPTS, MEQNS, NDMDS, X, Y, Z, 

IIA) 

Generate the interpolation information array for each modal input 
zone and save on MISFILE. 


(41) SUBROUTINE IIARDR (IIAP, IIA, LIIA) 
Read the interpolation information arrays. 


(42) OVERLAY (RHOIV, 1, 1) PROGRAM INPREP 
Process all user input. 


(43) SUBROUTINE KRNFCT 

Calculate the desired kernel function as specified by KIND, e.g 
GXYY, NON-SING, SINGULAR. 


(44) SUBROUTINE LETERM (GXIETA) 

Calculate the pressure term G (XI, ETA) for a leading edge control 
surface. 


(45) SUBROUTINE MATIO (LFN, MATRIX, MROWD, MROW, MCOL, LID, ID, 

IRR) , 

Read/write a two record set consisting of an ID record describing 
a matrix or array, and the matrix record. 


(46) SUBROUTINE MIINCK (NROWZ, X, Y, IPOS, ZONE, Z, IIA, NIPTS, 

1 PLNERR, MI ERR) 

Read modal input (from INPUT or MIFILE):, check legality, and 
write by input zone to INSFILE. Alternatively, read interpolation 
information arrays by input zone and write directly -to MISFILE. 

The information is used in MIPREP and RESPREP. 


(47) OVERLAY (PHOIV, 1, 2) PROGRAM MIPREP 

Allocate working area for the modal preparation routines which 
will generate information on MISFILE to be used in RESPREP. 


(48) SUBROUTINE NAMBLD (NMFILE) 

Convert NMFILE = NM to NMFILE = TAPENM 


(49) SUBROUTINE PCMDWM (NDWP, NPTRM, NOCS , NPMDS, C, CS , AS , 

W, WRI) 

Print the C-matrices associated with regular (main surface) 
assumed pressure terms and any control surface pressure term. 
Print the kinematic downwash matrix, and the residual downwash 



matrix (the kinematic downwash matrik with any control surface 
singularities removed) ... 

(50) SUBROUTINE PLATEO (XO, YO, NOPTS , ZO, NROWZ, NCOL1 , NCOLS, SA, 

1 INBD, DZ1 , DZ2) 

Given the spline coefficients for a set of functions as determined 
in routine PLATEI, and the associated input points, calculate 
the values of the functions (and optionally the derivatives) 
at a set of output points. 

(51) SUBROUTINE PCMSPT (NDWP, NPMDS, CMSPT) 

Print the coefficients of the regular (main surface) assumed 
pressure terms . 

(52) SUBROUTINE PLATEA (X, Y, INDS, SK, N, M, A, IRR) 

Form the coefficient matrix for system of equation, SK = smoothing 
constant (ratio of plate stiffness to input point spring 
. stiffness) . _ 

* 1 X(1) Y (1 ) Where A (I, J = R**2 LN(R**2, IAJ 

* - 

A (I , J) * - - = 0, I=J INDS = 0 

* - 

* I X (N) Y (N) = SK (1 ) , I=J INDS =1 

* 


= SK (I) , I=J INDS=2 


1 

1 

* I 

* 0 

0 

0 

R= (X(I) -X (J) **2 

X(1) 

— ~ X(N) 

* 0 

0 

0 

* (Y (I) -Y (J) ) **2 

Y (1) 

Y (N) 

* 0 

0 

0 

N=No. pts., M=N+3= No. eqs. 

(53) 

SUBROUTINE 

1 PLATE I 

(XI, YI, MIPTS, 

ZI, NROWZ, MCOL1, MCOLN, 


1 MCOLS , SA, INDS, SK) 

Perform a bivariate interpolation using as the interpolating 
function the small deflection equation of an infinite pinned 
plate. 

Reference: Robert L. Harder, Robert N. Desmarais; Interpolation 
Using Surface Splines, J. Aircraft, Vol 9 No. 2, Feb. 1972 
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(54) SUBROUTINE PLATET (XU, YV, NIND, XBAR, YBAR, COST,, SINT., 

rgu , rgv) V. ■ 

Perform a transformation of coordinates frcir (X,Y) to (U,V) if 
NIND>0, or from (U,V) to (X,Y) 


*0* * 1/RGU 0 * * COST SINT* *X-XBAR* •" , 

♦ ’ * =* . * * * * ■*_ ‘ 

*V* * 0 1 /RGV * *-SINT COST* *Y-YEAR» 


Where COST = COS (THETA) , SINT = SIN (THETA) , THETA is that angle 
such that PUV = SUM[ U (I) *V (I) ]N := 0, RGU, RGV = radii of qyratiea 
in (U,V) , XBAR, YBAR = C.G. location in (X,Y) - note UEAR^VBAR 
= 0,0 •' ' 


(55) SUBROUTINE’ PLNINCK (PLNERR) 

Read the leading edge and trailing edge definition, and check 
for compatibility. 


(56) SUBROUTINE PRSINCK (YPRC, XPPT, PLNERR, PRSERR) 

I . - • . ' : 

Read (or use default) pressure output chord and point values, 

check for legality and save on INSFILE. 


(57) SUBROUTINE PRSPREP (NPRC, NPPRC, NPPT, NSPT, . NCPT, NPTRM, 

1 MOCS, NPMDS, YPRC, XPRT, MSPTRM, FETA, CCR) 

Prepare information, independent of k-value and Mach no. , which 
will be used during calculation of unsteady delta pressures. 

The information is written to RESFILE. 

(58) OVERLAY (RHOIV, 1 ,4) PROGRAM REWRTC 


Perform all read/write activities associated with the CMFILE 
library of C-matrices and its index. 


(59) OVERLAY (RHOIV, 1,3) PROGRAM RESFREP 

Allocate working area for the various preparation routines which 
will in turn generate information which is placed on RESFILE 
for subsequent use in solving for the coefficients of the assumed 
main surface pressure terms, and in calculation unsteady 
pressures, sectional forces, or total generalzied force results. 



(60) OVERLAY (RHOIV, 1,7) PROGRAM RESULTS 

Allocate working area for the result routines which will in turn 
generate the requested results for printed or user file output. 

(61) OVERLAY (RHOIV, 1,0) PROGRAM RHOIV 

Calculate the unsteady aerodynamic loadings on a lifting surface' 
with leading and/or trailing edge sealed gap controls undergoing 
harmonic motions in subsonic compressible flow. 

(62) SUBROUTINE RTHETA (PRESS, AMPPHAS, NPPT, NMDS) ;; 

Convert pressure at a point from (X+IY) form to (R, THETA). 

(63) SUBROUTINE SCAMP 4 (X,Y,N,NDA,NDB,DA,DB,C,S,ML) 

Given a set of N points whose absicissae form a strictly monotone 
sequence, and given a first or second derivative at X (1) and 
a first or second derivative at X (N) , to find the smoothest 
possible curve passing rigorously through the given points,, 
satisfying the specified boundary conditions, and possessing 
continuous first and second derivatives. The criterion of i 
smoothness is the minimization of the integral of the square 
of the second derivative, and the curve found is accordingly 
a chain of cubics, i.e., a separate cubic on each interval X (I) , 
X(l+1) . 


(64) SUBROUTINE SERIES 1 

Routine SERIES 1 performs the integration of 

EXP (I*LAMEDA) / (LAMBDA * *2 + (K*YO) ** 2) ** 1 . 5 

In LAMBDA from - infinity to - K*F, where K*F = 2*PI. A change 
of variables is made, and a series expansion of the denominator 
is written. Using partial integration a recursion formula is 
developed which converges to within required accuracy in ten 
iterations or less. 


(65) SUBROUTINE SERIES 2 

Routine SERIES2 performs the integration of 


3/2 

(EXP (I*LAMBDA) -1 -I*LAMBDA + . 5*LAMEDA**2) / (LAMBDA** 2 + (K*Y0) **2) 


In LAM EDA from -K*F to K*H, when K*H > -K*F where K*F = 2*PI 
A series expansion of the exponent is written. The first three 
terms of the series cancel with the remaining three terms in 
the numerator. Performing partial integration twice# a recursion 
formula for the integral of the remaining terms is developed 
which converges to within the required accuracy within 50 j-.-.Vr. . 

iterations or less (K*H<1 5) . 


(66) SUBROUTINE SERIES 3 

Routine SERIES3 performs the integration of 

EXP (I*LAMBDA) / (LAMBDA* *2 + (K*Y0) **2**1 .5 

In LAMBDA from - infinity to K*H when K*H < -K*F, K*F = 2*PI. 
SERIES3 is the same as SERIES 1 except the sine and cosine 
integrals of the upper integration limit cannot be precalculated^ 

(67) SUBROUTINE SGFINCK (YSGFC, PLNERR, SGFERR) 

Read (or use default values for) sectional force output chords 
for legality and save on INSFILE. 

(68) SUBROUTINE SGFPREP (NSGFC# MPCHD, NZONES, NDMDS, NPMDS, NSPT, 

1 NCPT , NPTRM, YSGFC, NIPTSZ, XIPT, YIPT# 

2 CQWT, QTYPE, QSMSPT, Z, WZ, FETA, IIA, IIAP, 

3 CCR, MIPSGF) 

Prepare information, independent of k-value and Mach no.# which 
will be used for calculation of sectional generalized unsteady 
aerodynamic coefficients (sectional forces) . The information 
is written to RESFILE. - 


(69) SUBROUTINE SGRID (NIR, ETAQ# IRTYPE) 

Routine SGRID defines the spanwise integration scheme used to 
integrate the spanwise integrand of the downwash integral equation 
for all downwash points on a specified downwash chord. 


(70) SUBROUTINE SI Cl 

Computes the sine and cosine integral. 


(71) OVERLAY (RHOIV, • 1 #6) PROGRAM SOLN 



Allocate working area for the 
output routines. • • 


solution routine and intermediate 


(72) SUBROUTINE SOLVE (NDWP, NPTRM, NQCS, NPMDS, MMDS, 

1 C, SCR, IPR, GS, AS, W, WRI) 

Solve the coirplex linear system of equations for the coefficients 
of the assumed regular (main surface) pressure terms and save 
on COFFILE. 


(73) SUBROUTINE SPFCTE ' 

Calculate NSPT spanwise pressure terns 

KIND = GXYY, calculate F (ETAE) and FP (ETAB) - DF (ETAE) /DETA 
= NON-SING, or SINGULAR, calculate F (ETAB) 

CPT =.F. - main surface terms, A.T. - control surface terms 

(74) SUBROUTINE SPNINT 

Evaluate the spanwise integrand at ETA, and increment the C- 
raatrix terms fay the weighted results. 

(75) SUBROUTINE TETERM (GR, GI) 

Calculate the pressure term G (XI, ETA) = CMPLX(GR,GI) for a 
trailing edge control. 

(76) SUBROUTINE VPINCK (VPERR) 

Read velocity profiles, check for legality, form a cubic spline 
for each profile, and save on INSFILE. 

(77) SUBROUTINE ZEROCOL (M, NF, NL, Z, NROWZ, INDD, DZ1 , DZ2) 
Initialize columns NF-NL to zero for M rows. 


(78) OVERLAY (RHOIV, 0,0) PROGRAM ZRHOIV 
Call overlay RHOIV (1,0) 
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APPENDIX A - NUMERICAL INTEGRATION TECHNIQUES 

The basic work performed by the RHGIV program is numerical 
integration. In particular this includes the surface integration 
of the pressure kernel expression of the downwash integral 
equation (eqn 2.3-1) , or the integration along a chord of unsteady 
pressure times modal deflection for sectional generalized forces 
{eqn 2.8-4 ) , or the surface integration of unsteady pressure 
times modal deflection for total generalized forces (eqn 2.8- 
9) . The numerical integration is accomplished using Gaussian 
quadrature of the general form, 

e n . - 

/ w(t)f(t)dt, = (A. 1-1) 

a 


where f (t) = the function to be integrated with the associated term, 
w(t) = a known positive function of t, the weight function, 
and - the quadrature weights corresponding to, ~ 

a^ = the quadrature abscissae. 

The abscissae are roots of an nth degree polynomial which is 
orthogonal with respect to the weight function w { t) on (a,B) . 

A discussion of the existence and determination of the set of 
a^ and corresponding Fh for which eqn A.1-T will be exact, 
provided f (t) is of degree < 2n-1, is given in reference 3, 

Chapter VII. The properties of a- and include, 

a^ real and distinct 

> 0, i=1,...,n 

a < a^ < 8, i=1 ,. . . ,n 

n B 

E H i = / w(t)dt 
1 a 

The usual approach is to determine a. and H. for a known interval, 
e.g. (0,1), and perform a linear transformation of variables 
from, for example, x on (a,b) to t on (0,1) . 

The weight functions used within RHOIV are: 

w(t) = 1 (for which the abscissae are roots of an nth 
degree Legendre polynomial on (a, 8) — . 
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referred to as Legendre quadrature.) : 

/t-a (referred to as square root quadrature} 
l//t -a (referred to as inverse square quadrature! 
ln(t-d) (refered to as log quadrature! 


For integration performed in evaluation of the downwash integral 
equation, the inbegrand and its first derivative are- well behaved 
over the majority of the surface, thus Gauss-Legendre quadrature 
is most prevalently used. However, the first chordwise pressure 
term, cot (0/2) , for a main surface analysis has a square root 
singularity at the planform leading edge, and the remaining 
chordwise pressure terms for both main surface and control surface 
analyses have a J «qaar e root (or similar order) discontinuity 
in the first derivative at the planform leading and trailing 
edges. Additionally the spanwise terms for both analyses have 
a square root (or similar order) discontinuity in the first v 
derivative at the planform tips. The kernel expression contains 
a logarithmic singularity at a downwash chord (y) in spanwise 
integration, and a logarithmic singularity exists in the chordwise 
pressure form at a physical control surface hinge line. Thus 
Log quadrature is used in regions which include the above 
singularities. Note that because of the predetermination of 
the weights and abscissae on an interval (0,1) , the integration 
of a logarithmic singularity involves both Log and Legendre 
quadratures. 

For integration performed in calculation of sectional and total 
generalized forces the kernel function singularities are not 
involved, however the singularities associated with the pressure 
terms remain and are handled in the same manner. 

The form of the quadratures used is. 


w(t) = 
w(t) = 
w(t) = 



Gauss-Leqendre Quadrature 


b 1 n 

/f(x)dx = (b-a)/f(t)dt = (b-a)EH^ -f (x^ ) 
a 0 1 


(A. 1-2) 


where n = 4 or 8 

the n point Gauss -Legendre quadrature weights 

x - = (b-a)t. + a, t. = the n point Gauss-Legendre 

quadrature abscissae (note 
that H. and x. are symmetric 
about ^b+a) /2 1 

S quare Root Quadratur e 

b 1 1 

f(x)d X = (b-a)//(b-a)t f(t)dt = (/tT^T ) 3 f/t f(t)dt (A. 1-3) 
a 0 o 

= (/E^a ) 3 zHSif(x,. ) 

1 1 

where n = 5 or 10 

Hs^ = the n point Square Root Quadrature weights 

X. = (b-a)t; + a, t. = the n point Square Root 

1 1 Quadrature abscissae 


f (x) = g (x) / /x-a, where g (x) has the characteristic of 

/x-a 

Inverse Squ a re Root Quadrature 

/f(x)//x-a dx = (b-a)/f (t)//(b-a)t dt = /b-a /f(t)//t dt (A. 1-4) 
a 0 0 

n 

= /b-a EHis.-f(x-) 
l 1 i 


where n = 5 or 10 
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His^= the n point Inverse Square Root Quadrature weights 

x. = (b-a)t. .+' a, t. = the n point Inverse Square Root 
1 1 1 Quadrature abscissae 

f (x) = g (x) /Fa/ where g(x) has a sq uare root singularity 

of the form 1/ /x-a 


Log ( Plus L egendre) Q uadratur e 

b 1 

/In] x-a | f(x)dx = (b-a)/1n| ( b-a ) t | f(t)dt 


(A. 1-5) 


1 • : 1 \ 

= (b-a) 1 n | b-a|/f(t)dt + (b-a)/lnjt| f(t)dt 
0 0 
n n . . 

= (b-a)1n|b-ajiH.f(x.) - (b-a)EHl -f (x .) 

iv i — J . 

where n = 4 or 8 

H. and x. are as defined in C.1-2 
1 i 

Hlj= the n point Log quadrature weights associated with 
w (t) = 1 n (1/t) • 

x. = (b-a) t . + a, t. = the n point Log quadrature 
3 3 3 abscissae 


f (x) = g (x) / In | x-a |, where g(x) has a logarithmic 

singularity of the form 1n|x-a| 
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APPENDIX B - MODAL INTERPOLATION 


If input points and modal displacements are provided by the user 
to RHOIV, the surface spline method of Harder and Desirarais, 
Reference 5, is used to allow interpolation for displacements 
and streamwise slopes at downwash points, and displacements at 
quadrature points used for sectional and total generalized force 
integration. 

If control surfaces exist, the lifting surface planform is divided 
into modal input zones, as illustrated in figure 8, page 155, each 
of which must have a sufficient number of input points to define 
the motion within that region. 

The user has two options when describing lifting surfaces with 
controls: 

(1) Provide a descrete set of input points with associated modal 
displacements for each modal input zone, 

(2) Provide a total set of points with associated modal 
displacements which includes points in all zones. 

In the first case, the input points for an input zone need not 
lie within the boundaries of that zone. In the second case RHOIV 
will assign an input zone to each point based on the boundaries 
shown in figure 8. In either case a minimum cf three (3) input 
points must be associated with each input zone; the interpolation 
procedure described below, is applied per input zone . 

The control surface pressure modes are a function of control 
rotation, 9 (n) , or the change in streamwise slope across a control 
hinge. Unless the user specifically provides this information, 
RHOIV will determine 9 by calculating the slope at the hinge 
using the interpolating functions for the two medal input zones 
on either side of the hinge. 

With a reasonably even distribution of points in a zone, the 
surface spline approach provides good results for displacements, 
and reasonable results for slopes at points interio r to the input 
point set. At the extremeties of input point regions a fair 
amount of curvature may exist, introducing sometimes large 
"errors” in the slopes. This is particularly significant since 
the leading edge and any control hinge tend to lie at the 
extremeties of input point sets and thus errors in boundary 
condition at the leading edge, and control rotation 9, may be 
introduced. 
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For any particular planform configuration, the user should 
initially examine the kinematic downwash and control rotations 
to assure the sufficiency of his input point distribution. If 
the user is unable to cause reasonable control rotation values 
to be calculated by moving input points or including more input 
points the option to input control rotation information should 
be used. 


IF REASONABLE SLOPES AT DOWNWASH P OIN TS. AND CONTROL ROT ATI ON 
VALUE S ARE NOT CALCULATED THE PROGRAM RESULTS WILL HAVE NO 

significance " 

If the surface spline approach is used, the modes are approximated 
'••‘^yby , (Reference 5), 


N 


Z(C/H) - j a^R^ • In (R^ ) + a^ + ^ + + 


R i 


= (C-x i ) 2 + (n-y i ) 
= N input points 


under the constraints 
Z(x i/ y i ) = Z l (x i /Y i >/ 1 = 1,N 

where Z_(x. ,y.) = input values at (x. ,y . ) , and 

^ j • ;X X 

j ■ • . ■- 

N N ' N 

£a. =0, Za.x. = 0, la.y . =0 
I r ill r ri ■ 

which expressed in matrix form for n input modes. 


1 x 


1 y l 


A ij 


1 . 

x r 

y r 


.rr 

,X nI 


• • . • 


*N y N 


r N 


c { a.} L . . .{a} n l = [:{Z I (x i ,y i ) > 1 ... {Z I (x i ,y i ) 




•( 


0, i=j 

Rfj -ln(R?j) , i^j 
(x i -x j ) 2 + (y i ~Yj ) 
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is a linear system of equations which may be solved for the 
interpolating function coefficients. 

Then for any output point, (£ ,n ) , 


[z (£ #tv) i f • • • g z (£ #n) n J ~ [s^j] J ^ 


= R?*ln (R?) , i 


B i = 2(C-x i ) [In (R?) +1] , , 


i=l,N 

where Z {% ,n) has the properties 

z < x i'Yi> = Z I ^ x i'Yi) f i = 1/N 

9 2 Z (5/h) 3 2 z (£»n) 9 2 Z (g , n) -» 0 as the distance to the 

3X 2 ~ ' 9x9y dP input point -*> » 

Z(C,h) is analytic everywhere except at (x^,y^) 

9Z(£,n) 9Z(C,n) exist at (x.,y.) 

53E ' *9y 1 1 


If N = 3 and the points are non colinear, a simple plane is 
defined for Z {£ ,n ) . 

Note that a simple bending mode (1st, 2nd, etc.) is n ot adequately 
defined by input points lying in two lines because of the tendency 
of the interpolation function used to force curvature to zero 
in both coordinates. 
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The change in slope at a control hinge for a spanwise station 
y is determined by calculating the slope at the hinge using first 
the interpolation information for the surface to which the control 
is related and then using the interpolation information for 
control surface itself. Rather than calculate the slope 
analytically, which, in the case of a surface spline 
representation, may introduce control rotations where none exist, 
a numerical procedure is used. 

At a spanwise station y, three equally spaced points (including 
the hinge) are selected over a sufficiently small region, S, 
e.g. 26 = (x 3 -x t ) = . 005S, so that it can be assumed that no 
curvature of the physical planform exists. If the displacements 
at the points are determined, a quadratic in (x-x t ) may be 
written, 

Z (x) ♦ A 1 +B(x-x l ) + C(x-x l ) 2 , or 
3Z(x)/3x = B + 2C(x-x t ) 2 

Solving the system of three equations and neglecting C, 

3Z(x)/3x = B = [ 4Z <x 2 ) - 3Z(x x ) - Z (x 3 ) J/2 6 

Note that if the slope is a constant over (x r , x 3 ) initially, 

c. * [ (Z(jc 3 )-Z(x,))/ 6 - (Z(x 2 )-Z( Xl ) )/6 ]/6 

This procedure should reduce the introduction of extraneous 
control rotations in non control rotation modes. 
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